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1 . OBJECTIVES AND SCOPE 


The objective of this project Is to develop a realistic 
biotechnology assessment tool. Its main benefit is seen as 
providing ECDT and industrial chemical manufacturers with a con- 
sistent basis for examining proposed bioroutes to bulk chemicals. 
The basis of this tool would be currently available technology but 
would allow asking 'what if' types of questions and thereby ascer- 
tain the total potential impact of new technology were it to be 
developed in either energy conservation or biochemistry (recombi- 
nant DNA related types of technology) . During the course of these 
evaluations* the energy requirement of bioprocesses, including 
upstream processing (raw material preparation), sterilization, 
media preparation, fermentor (reactor) requirements, and down- 
stream processing, which includes cell-liquid separations, solute- 
liquid separations, dryers and packaging and various interactions 
which can occur when 'what if types of questions are asked, will 
become more clearly understood and documented. This will be of 
value to ECUT in deciding which areas of bioprocess technology 
would most usefully profit from government-sponsored research and 
yet constitute too high a risk for industry to sponsor alone. 

DOE/ECUT is concerned with bulk commodity applications of 
bioprocesses. Considerably less emphasis should be placed on 
bioroutes to new pharmaceuticals and small-scale specialty 
chemicals whose high value-added market potential more than justi- 
fies any unit cost of energy. The project directs computer 
analysis at large scale processing, where even small but signifi- 
cant savings of energy on a unit basis have a broad impact on the 
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national energy budget. 

What is being offered is an accurate tool for consistent 
analysis of bioprocesses across a spectrum of potential appli- 
cations. This assessment can yield precise information about 
where DOE should encourage new research to provide the nation with 
large scale renewable raw material resources and energy- 
conservative technology. It opens up a gauge of prospective 
success for any high technology group wishing to evaluate the 
results of genetic engineering research. 

Industry is too constrained by short-range concerns to under- 
take the long-range research which will bring the potential of 
bioprocessing of renewable resources to fruition. Witness the 
recent abandonment of major synfuels projects by industry when the 
price of petroleum stabilized. Nevertheless, it is in the 
national interest that replacement technology be available 10 or 
20 years from now. Non-renewable raw materials will dry up inevi- 
tably. The potential of bioprocessing to make even non-renewable 
materials more attractive than conventional processing in the long 
run is not discounted. 

The System for Biotechnology Assessment is a small investment 
in the future. The cost is greatly transcended by its promise of 
rapidly screening opportunities and goals for DOE/ECUT. This tool 
would also serve as a valuable training and screening facility 
both in academia and industry. 

The specific objective is to convert the existing computer 
program known as DOE-ASPEN to function for bioprocesses. ASPEN 
was developed successfully to allow assessment of synfuel 
projects. The changes would add the unit operations peculiar to 
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bioprocesses and bioseparations, such as fermentation. The data 
base would be extended to include bacterial and cell properties. 
The costing routines would be updated and extended. To validate 
the system the already documented acetone-butanol process would be 
examined, as well as selected other bioprocesses of interest to 
ECUT/DOE. It is anticipated that towards the middle of the 
project a consortium of industrial companies would be formed to 
continue and refine the development of this tool and to facilitate 
its transfer to industry. This group should include represen- 
tatives of both large and small firms. 

A second, and unfunded, phase of the project is contemplated 
to begin in the later stages of the work. This would analyze the 
potential of developing a single technology assessment tool that 
can be executed on a user-friendly microcomputer such as an IBM 
Personal Computer. The downsizing of the large model would have 
to be done with some attention to the level of confidence 
attainable with restricted models. 
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2 . INTRODU CT ION 


The chemical and petroleum processing industry is faced with 
increasing shortage of raw material and increased costs of produc- 
tion over the next two decades. As petroleum supplies dwindle it 
will be more and more dependent on renewable resources. The new 
processes, based on biocatalyzed conversions, hold a high poten- 
tial of providing an alternative source of raw material and energy 
for the n.S. chemical industry. Although these processes are not 
currently competitive with other conventional sources, further 
research into new catalysts, new routes and improved engineering 
(reactors and separation processes) for production are increasing 
their potential. 

The question naturally arises as to what processes hold the 
greatest potential and what areas are most critical in terms of 
research needs. One would like to concentrate available resour- 
ces on those concepts which hold the greatest opportunity for 
improved efficiency and production of alternative fuels for use in 
energy conversion systems. 

At present no form of computer-based bioprocess analysis 
system exists. There is a need for such a system for the follow- 
ing reasons: 

(1) Bioprocesses are currently perceived as a major opportunity 
in the chemical and process industries 

(2) Different groups of investigators arrive at different assess- 
ments of whether a particular process is economic or not 

(3) Such assessments of bioprocesses are of necessity limited and 
use different starting assumptions 

( 4 ) DOE/ECUT is faced with information of this type and has few 
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methods of assessing its worth. It then has to decide which, if 
any, biotransformations are worth supporting. In particular, it 
must attempt to identify the weak links within a process, 
strengthening of which would allow the process to succeed. 

A TOOL FOR TECHNOLOGY ASSESSMENT is needed rather than a 
plant design tool. However, the degree of detail needed in 
process analysis for a realistic estimate of the worth of a bio- 
process is considerable. Part of the problem with the type of 
analysis available from consulting firms is that it inevitably can 
only provide answers to within 30-50% accuracy. This degree of 
doubt is rarely, if ever, made explicit. It is not surprising 
that industry, knowing these limitations, will be slow to move 
even on an optimistic forecast. However, recognition of, for 
example, the advantages of a 2-5% improvement in yield could in 
many cases be sufficient to persuade a bulk producer to switch 
to a bioroute for a desired product. The producer therefore 
needs a sufficiently accurate analysis to enable him to make that 
decision. 

It is anticipated that the models developed will have high 
comparative accuracy and a high degree of uniformity. They will 
allow, 1) direct comparison with conventional routes to the same 
chemicals, and 2) ranking of a range of proposed bioprocesses in 
order of likely economic return. This will permit ECUT to iden- 
tify research needs and targets more closely . Absolute accuracy 
in such models is less necessary and will depend largely on the 
detail of the simulation. 

(5) The system will allow individual companies to ask 'should 
chemical A be made by a biological route?'. The answer will 
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depend on who is asking the question, when, and for which loca- 
tion. Each company has its own specific special circumstances 
that it is usually unwilling to discuss publicly. These may 
include special feedstock agreements at favorable prices, conve- 
nient access to markets or raw material sources: a company may 
already have byproducts available from one plant that can be 
utilized in another, or special Interest rates and risk factors. 
Furthermore each variable will be subject to change as the company 
grows or diversifies and as the external environment changes, e.g. 
oil prices, and as management attitude to perceived risk changes. 

None of the above factors can be influenced directly by ECUT. 
They can be turned to advantage by providing a technology assess- 
ment tool that allows the question to be investigated directly by 
the end user, with explicit regard to his particular circum- 
stances, but without either the circumstances or the answer being 
disclosed to outside parties. 

(6) Much of the technical thrust for the development of biotech- 
nology has come from small companies consisting of outstanding 
molecular biologists, etc. Their main Interests do not and, at 
least for now, cannot be in the development of sophisticated 
economic and technology assessment tools. The proposed 
technology assessment tool would permit them to identify and 
exploit relevant research and commercial goals. 

Despite the requirements of a mainframe computer, the program 
would be directly available through timesharing systems over the 
telephone lines, and hence directly available to smaller 
companies. 

(7) The development of the technology assessment programs would 
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constitute valuable training for the graduate students concerned. 
After graduation they would be well placed to transfer this 
training to an appropriate biotechnology-based industry. The 
availabilty of the programs would also offer many institutions the 
opportunity to enhance their educational and training activities 
in this most valuable component of the nation's future economic 
well-being. 
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3. SUMMARY OF TECHNICAL PROGRESS 


The public version of ASPEN available through the Argonne 
National Energy software center was installed .in the Chemical 
Engineering Department VAX 11/750 computer. To examine the idea 
of BIOASPEN, a test example (the manufacture of acetone, butanol 
and ethanol through a biological route) was chosen for simulation. 
Previous reports (1, 2) on the BIOASPEN project have shown the 
limitations of ASPEN in modeling this process. 

To overcome some of the difficulties, modules were written 
for the acid and enzyme hydrolyzers, the fermentor, and a steri- 
lizer. Information required for these modules was obtained from 
the literature whenever available. Additional support modules 
necessary for Interfacing with ASPEN were also written. Some of 
the ASPEN subroutines were themselves altered in order to ensure 
the correct running of the simulation program. After testing of 
these additions and changes was completed, the Acetone-Butanol- 
Ethanol (ABE) Process was simulated. 

A new release of ASPEN (which contained the Economic Sub- 
system) was obtained and installed. This subsection was tested, 
and numerous changes were made in the FORTRAN code. Capital 
investment and operating cost studies were performed on the ABE 
Process. Some alternatives in certain steps of the ABE simulation 
were investigated in order to elucidate their effects on the 
overall economics of the process. 



^ DETAILED DESCRIPTION OF TECHNICAL PROGRESS 


This section describes in detail the progress made in the 
various aspects of the project to date. 

4.1 Installation of ASPEN 

The public version of ASPEN was installed in the Chemical 
Engineering VAX computer. The test examples provided by the 
Morgantown Energy Center were run through the program. UOSD, the 
test problem for the solids handling modules, could not be run 
because of the lack of sufficient documentation. Work was 
therefore started on the development of a documentation aid for 
ASPEN programmers, which would give information on the hierarchy 
of subroutines used in any process module. With the installation 
of the later version of ASPEN (version 6 ) , this work was 
discontinued as an on-line documentation aid had been provided 
with the package. 

4.2 Incorporation of New Models into ASPEN 

ASPEN, as is. was considered inadequate to model various 
aspects of the ABE Process (1. 2). It was considered essential to 
model the acid and enzyme hydrolyzers. and the fermenter. A 
sterilizer was also incorporated into ASPEN. While the actual 
FORTRAN code for each of the modules may be found in the Appendix, 
this subsection provides information on the actual models used in 
each of these units. 

(i) Acid Hydrolysis: 

The kinetic model in general use (3. 4. 5. 6 etc.) for the 
cellulose hydrolysis seems to be the model proposed by Saeman (7): 
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crystalline cellulose > glucose > degradation products 

amorphous cellulose > glucose > degradation products 

Chambers (8) proposed a similar model for hemicellulose 
hydrolysis: 

hemicellulose — -> xylose > degradation products 

The rate constants k^, follow the Arrhenius rate form: 

k^ ■ A^exp(-Ej^/RT) 

with the pre-exponential factor being a function of acid con- 
centration as follows: 

*1 ■ *ol=" 

The kinetic constants used in this work were reported by 
Bhandari et al (6). They make the assumption that the amorphous 
cellulose* which is the easiest fraction of cellulose to hydro- 
lyze, is attacked and converted instantaneously to glucose. 
Therefore, an initial yield of glucose was assumed equal to the 
amount produced from the amorphous cellulose. They made a number 
of experimental runs, and produced best-fit equations for cellu- 
lose and hemicellulose hydrolysis reaction constants based on corn 
stover as substrate: 

kj - 2.7*10^^C^*'^^exp(-45,300/RT) 
kj - 2.01*10^^C^**^exp(-32,800/RT) 
kg - (7.84 - 3.68/C)*10'‘”exp(-41.000/RT) 
k^ » (4.80 - 1.95/C)*10^^exp(-32,000/RT) 
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They also present a table comparing their kinetic constants 
with those available in the literature for other substrates. Some 
of this data along with kinetic constants presented by Kwarteng 
(5) for Aspen wood are shown in Table 4.1. and are Incorporated 
into the user model. 


(ii) Enzyme Hydrolysis: 

The literature provides a number of models that may be used 
for the enzymatic hydrolysis of cellulose (9, 10, 11, 12 etc.). 

Most of them however, are applicable to restricted systems (speci- 
fied substrates and microorganisms), or require enormous numbers 
of parameters to be specified. In order to keep the model simple 
at this time, it was decided to use rate expressions presented by 
Pietersen et al (12). This model does make a large number of 
assumptions, but also contains elements of the common observables 
of enzymatic hydrolysis of cellulosic substrates: 

(1) Cellulose consists of a crystalline and an amorphous 
region. 

(2) Michaelis-Menten kinetics may be used. 

(3) Reactions are subject to product inhibition. 

(4) The organism grows on glucose, with a Monod growth rate. 

The reaction scheme is as follows: 

crystalline cellulose > amorphous cellulose > cellobiose 


\ X 


glucose 

The different rate expressions (in g/liter hr) are: 
i) crystalline cellulose to amorphous cellulose 


rate - 0.959 f [Cl [E ^l 1 g/liter hr 

1.2.14 + [Cl +^46.5[AlJ 
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Table 4,1 Comparison of Rate Parameters for Acid Hydrolysis 
(Different Substrate Materials) 



Preexp. 

Factor 

Substrate 

Material 

(min 

•icr* 

) 

•10'’^ 

Douglas Fir 

17.3 

2.38 

Kraft Paper 

280 

4.9 

Solka Floe 

12.2 

3.79 

Oak Sawdust 

4.4 

0.028 

Corn Stover 

27.1 

2.01 

Aspen 

1.45D-3 

3.84D-5 


Acid Activation 


Concentration 

Energy 



(kLn/9mo!f* 

1.34 

1.02 

42.9 

32.87 

1.78 

0.55 

45.1 

32.8 

1.16 

0.69 

45.2 

32.8 

1.00 

1.80 

42.9 

30.0 

2.74 

1.86 

45.3 

32.8 

1.16 

0.57 

33.7 

21.0 
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ii) amorphous cellulose to cellobiose 


rate - 31.3 f [A]^*^ [E^l 1 g/ liter hr 

L 30.4 + [A] + 9.74[BlJ 

lii) amorphous cellulose to glucose 

rate » 4.12 r [A1 [E-] 1 g/liter hr 

[44.1 + [A1 +''2.25[GlJ 

iv) cellobiose to glucose 

rate « 3.58 f [Bl LE- 1 ^Ig/liter hr 

L0.265 + [Bl +‘^0.376[G1J 

Here, concentrations are in g/liter and the enzyme concen- 
trations [E^l and [Ejjl are in FPU/ml (filter paper units of 
cellulase activity) and BD/ml (cellobiase units of P-glucosidase 
activity) . The microorganism used in their study was T. Reesei . 
QM 9414, which is the same as that used in the present model. 

(iii) Acetone-Butanol-Ethanol Fermentation: 

For an anaerobic fermentation, a balance equation may be 
written: 

CH^Oi + aNHg -> ycCHpO^Nq + zCH^O^N^ + bH^O + cCOj 

The problem however, is in determining the yield coefficients 
(a,b,c,y^, z) , and the kinetics of the reaction system. 

The reaction kinetics and the yield coefficients depend on 
the biochemical pathways available to the system. For butyric 
acid bacteria, these pathways for glucose fermentations are 
numerous (see Figure 4.1) and complex, involving a large number of 
intermediate/final compounds. Moreover, the kinetics of each step 
are not known. An analysis of the pathways however, yields impor- 
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Figure 4.1 Biochemical Pathway of Glucose 

Fermentation by Butyric Acid Bacteria 
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tant information on the yield coefficients. This topic has been 
addressed in detail by Papoutsakis and his coworkers (13. 14). and 
they present an overall equation for solvent production as 
follows : 

(2 + a) glucose — > + (2a - b - 2k) pyruvic acid + 

(b-e-f-g-h-j) AcCoA + (2a - 1,75 + 
d - 2f - 2g - i - 2j) NADHj + (b - C - d) FdHj 
+ (2a - 29,7 + e) ATP + (e + g - f - h) acetic 
acid + (f - g) butyric acid + cHj + (b + h + 
2k) COj + g butanol + (h - i) acetone + i iso- 
propanol j ethanol k acetoin 

To fully determine the above, the 11 constants a through k 
have to be determined. Some additional equations are valid for 
these fermentations, and they are listed below: 

2a - b - 2k » o’ 
b-e-f-g-h-j«o 
2a - 1.75 +d-2f-2g-i-2j-0 

b - c - d « 0 

The above equations are valid because AcCoA, NADH 2 , FdH 2 and 
pyruvic acid do not accumulate in the fermentation broth. These 
equations have been shown to include the Carbon balance and an 
available electron balance. Also, the derivation of the overall 
equation Includes the assumption that (i) the ATP yield * 

10.5 g/mol, and (ii) the fraction of carbon in biomass. T • 0.462 
* 0.023. The overall equation has been shown (13), based on data 
available in the literature, to be valid. 

A computer survey of the literature produced only one model 
for the ABE fermentation (15). The results reported were not 
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reproducible, possibly because of an error in some parameter. 
Therefore, in the absence of any kinetic model for the fermen- 
tation, a yield-based model utilizing the overall equation was 
used. For the purposes of BIOASPEN at the present time, kinetics 
are not as important as the knowledge of correct yield values of 
the products. This is because the yield coefficients have an 
immediate effect on the economics of the process. Therefore, a 
Monod growth model of biomass growth and substrate (glucose) 
utilization was used, with the product formation being linked 
directly to growth. It should be noted that this model would, 
therefore not show the production of acids (which are intermediate 
products in the biochemical pathways) prior to solvent formation. 

(iv) Sterilization: 

Sterile operation in bioprocessing is critically important 
and may account for a fair proportion of operating costs. Heat 
sterilization, whether continuous or batch, is also energy inten- 
sive. Extensive discussions on sterilization may be found 
elsewhere (1«). A short description of the model equations for 
the kinetics of thermal death is presented here. 

Logarithmic death rates may be represented by: 

g « exp(-kt) 

For bacterial spores, the non-logarithmic death rates that are 
found may be described by a sequence of reactions, in which the 
formula becomes: 
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The rate constants k, and k^ may be represented by the Arrhe- 
nius rate forms: 

k - A exp[- ||] 

Deindorfer and Humphrey (17) describe equations for batch sterili- 
zation. The del factor V is defined as: 

Vtot -Vheat *^hold *Vcooi * R - » «*p[- ||]at 

o 

The time-temperature relationships for heating and cooling depends 
on the method of heating/cooling, and are shown in Table 4.2. 
Depending on the values of the various parameters (such as heat 
transfer coefficients, temperature of sterilization, flow rates of 
steam) the steam requirements may be worked out. 

For continuous sterilization, one %ay use the continuous 
injection method or the continuous plate exchanger. In both 
cases, extra units are required in addition to the steam, adding 
to the capital costs. The material balance may be written as: 

isp - *■* a? - “ 

where N - N/N^ 

o 

Pe ■ Peclet number ■ UL/D„ 

z 

* axial dispersion coefficient 
X - X/L 

Da « Reaction number * kL/U 
U » average velocity 
L « length of the unit 
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Table 4.2 Temperature-Time Profiles in Batch Sterilization 


Type of heat transfer 

Temperature-time profile 

Parameters 


Steam sparging 

hyperbolic 

hs 

b. * 

"""mToPC, 

M 

Electrical heating 

T - To(1.0 -1- at) 
Imar 

a- ’ 


' MTopC, 


Steam (heat 
exchanger) 

r-TH(l + 6e-) 

ciponeniatl 

i''* I. 
‘-mT * 

1 

I 

Coolant (heat 
exchanger 

T-r^(l + fie-) 

expooeotMl 




where h - enthalpy diflerenoes between steam at sparger temperature and raw 
medium temperature 
5 « steam mass flow rate 
M “ initial medium mass 
7^ > initial medium temperature 
q « rate of heat transfer, kcai per unit time 
U « overall heat-transfer coefficient, kcal/m^‘h’*C) 

A » heat-transfer area, m^ 

Tm “ temperature of heat source 
w ■ coolant mass flow rate 
c* - coolant specific heat 
7^ « coolant inlet temperature 
p « medium density 
C, “ medium heat capacity 
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The boundary conditions are: 


X-> 0 


X-*- 1 


+ Pe(l - N) - 

dN 
dX " 


The solution to the differential equation is: 


III » exg(Pe/2) 
where y - [l ♦ 45|] ' 

The Peclet numbers may be estimated for fluids flowing in 
circular pipes at various fluid dynamic conditions. The length 
requirements may therefore be calculated, or alternatively, with a 
given length one may calculate N 

4.3 Testing of the Economic Subsection in ASPEN 

The new release of ASPEN (version 6) contained the costing 
modules missing from the previous release. Sample problems were 
formulated and run on this new system, to verify the smooth 
running of this subsection. Major problems surfaced during these 
exercises. These were all corrected, and the changes in the 
FORTRAN code for each subroutine were listed in a separate sub- 
directory. The original subroutines were also retained so that a 
proper account was kept of each change. 



4.4 Simulation of the ABE Process 

A flowsheet was formulated based on a report by Chem Systems 
Inc. (18), and ASPEN was used to simulate this. Figure 4.2 shows 
the broad subsystems in the process. The ABE Process basically 


20 


enzyme 



21 


Figure 4.2 Overall A - B - E Process 




consists of fermentation of sugars to solvents. A biological 
route to the formation of sugars from the cellulosic components of 
wood is also incorporated, resulting in the six-part process 
described below. 

(i) Prehydrolysis 

Acid hydrolysis under carefully controlled conditions (0.5 
percent H^SO^, 374 deg. P, 185 psia. 12 secs residence time in a 
plug flow reactor) causes the conversion of hemicelluloses and 
amorphous celluloses - which are the cellulosic components of wood 
- to glucose and xylose. Some of the cellulose is converted to 
the degradation products, hydroxymethyl l-2-furfuraldehyde and 
furfural. Wood chips are slurried in a prehydrolysis slurry tank 
where they are slurried to 35 percent solids, and conveyed to a 
plug flow reactor. The product stream is immediately quenched, to 
reduce the degradation reactions. The solids are separated from 
this stream and sent to the enzyme hydrolysis section, while the 
soluble sugars are sent to the fermentation section. 

(ii) Enzyme Hydrolysis 

The remaining unconverted cellulose (crystalline cellulose) 
is hydrolyzed enzymatically to glucose in this stage. The enzymes 
are produced by a mutation of T. Reesei . RUT-C-30. Typical con- 
ditions required for this stage are 122 deg. F. a pH of 4.8 and 
hydrolysis time of 24 hours. Conversions of cellulose to glucose 
may be as high as 90 mol percent . with terminal glucose concen- 
trations of the order of 5 percent. 

The solids from pretreatment are slurried and fed to the 
enzyme hydrolyzers. which are cone-roof atmospheric tanks with 
heating coils. The product stream is filtered, with the filtrate 
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going on to fermentation. The solids, consisting mainly of lignin 
and unconverted cellulose are burnt as fuel. 

(iii) Enzyme Production 

A mutation of the T. Reesei fungus, RUT-C-30, produces endo- 
gluconase and p-glucosidase, the enzymes required for enzyme 
hydrolysis. Continuous fermentation at 86 deg. F and a pH of 4.8 
is done with ligno-cellulose as a carbon source and corn-steep 
liquor as nitrogen source. Other nutrient requirements are provi- 
ded for by adding organic salts. The outlet from the fermentor is 
centrifuged: the filtrate is jsent to the enzyme hydrolysers, while 
a portion of the solids is recycled to the fermentor and the 
remainder is recovered as a single cell protein by-product. 

(iv) Fermentation 

Fermentation of the sugars is done with C. acetobutylicum. 
Since the microorganism cannot tolerate a butanol concentration 
greater than 1-3 percent, there is a limit to the initial sugar 
concentration (6 percent) that may be used. Experiments however, 
have shown that the yield on this fermentation begins to decrease 
as the initial sugar concentration is increased beyond 3 percent. 
Accordingly, the Chem Systems report studied two cases: one 

maximizing yield, and the other with a 5 percent initial sugar 
concentration. Since the latter case was shown to be more econo- 
mical, it was chosen for this study. 

A typical yield for this case is 30.6 percent on 90 percent 
sugar utilization after 48 hours fermentation time. The overall 
yield is thus 27.54 percent, with a solvent ratio of 61.7 percent 
butanol, 31.8 percent acetone and 6.5 percent ethanol. Corn is 
used as a growth medium for C. acetobutylicum culture maintenance, 
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and ammoniuin sulfate as the nitrogen source. Calcium carbonate 
and superphosphate are also required in small amounts. 

The fermentation vessels are cone-roof tanks in which the 
sugar streams from prehydrolysis and enzyme hydrolysis, the 
nutrient media and the seed inoculum are added during the fill 
period. COj and Hj evolved during fermentation is sent to the COj 
recovery section, and the fermentation beer is sent to a beer well 
prior to purification. 

(v) Purification 

In this section, the product stream from fermentation is 
split into pure acetone, butanol and ethanol. A sequence of 
distillation columns is used for this purpose. Steam costs in 
this section are extremely high, mainly because the outlet stream 
from fermentation is very dilute. 

(vi) Carbon Dioxide Recovery 

Because the amount of carbon dioxide produced during fermen- 
tation is large (as much as 3.3 ft^ per lb of sugar consumed), it 
is economical to recover this as liquefied CO 2 for sale. 

The off-gas is passed through a water scrubber to remove 
impurities and then compressed. This compressed gas is then 
passed through activated carbon columns to remove other impurities 
and is then chilled. The chilled gases then pass to a stripper- 
condenser system, where the CO 2 is separated from the other gases. 

A simulation was carried out using the user modules as 
described in Section 4.2. Some sample input data files and output 
reports can be found in the Appendix. The separation of the 
products from the fermentation broth involves fairly complex 
distillation column calculations, including azeotropic 
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calculations . The flowsheet presented by Strobel and Bader (19) 
uses an improved separation sequence (see Figure 4.3) than that 
used in the conventional I950's scheme (here shown in Figure 4,4). 
In this scheme, the butanol-water azeotrope is separated earlier. 
A beer column concentrates the broth to near the azeotropic 
composition and feeds it to a decanter. The butanol-rich phase is 
fed to a second column, from which a pure butanol stream is 
withdrawn. The third column does an azeotropic separation, with a 
butanol-water mixture being recycled back to the decanter. The 
acetone and ethanol are separated further downstream. The 
improvements of this process over the conventional one are: (i) 
the major product, butanol, is removed earlier, resulting in lower 
flow rates downstream, and (ii) there is one less column required 
in the overall scheme. 

This flowsheet was therefore used for the present simulation. 
Rigorous calculations however, were not performed for any of the 
distillation columns; instead, a sequence of splitters and mixers 
with defined splits (based on the calculations of Strobel and 
Bader) were used. For the costing exercises, the sizing calcu- 
lations performed by the same authors were used, as they con- 
sidered a product stream of a similar flow rate as that produced 
by the present ABE simulation. 

Figures 4.5 through 4.7 show the process flowsheet. Wood 
chips enter the mill (Figure 4.5), where they are crushed. The 
crushed wood is slurried in a tank, and then introduced into the 
acid hydrolyzer at the appropriate conditions. The quench serves 
to arrest the degradation reactions. After further cooling , the 
soluble sugars are separated from the crystalline cellulose and 
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Figure 4.3 Improved Separation Sequence 




Figure 4.4 Conventional Separation Scheme 
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Figure 4.7 Fermentation 


lignin , passed through a carbon column for purification purposes, 
and sent on to the fermentation section. 

The unconverted celluloses and lignin from acid hydrolysis is 
mixed with water, acid and enzyme , and introduced to the first of 
a train of enzyme hydrolyzers (Figure 4 , 6 ). These tanks are 
maintained at the appropriate conditions, and after a total 
residence time of about 50 hours, the product stream is filtered 
in a rotary filter. The filtrate is sent on to the fermentation 
section after removing the enzyme if desired, and the solids are 
used as fuel. 

In the fermentation section (Figure 4.7) the two streams from 
the two hydrolysis steps are sterilized by steam injection, and 
then cooled down to the temperature required for fermentation. 
The off-gas generated from the fermentors are compressed to 
facilitate carbon dioxide removal. The liquid product stream is 
sent to the purification section (see Figure 4.3), and the 
recycle is fed back into the reflux mixer. A portion of this 
stream is also purged, so that there is no buildup of biomass. 

4.5 Costing of the ABE Process 

ASPEN uses cost correlations which were produced by reg- 
ression on cost data for a large number of equipment classes. 
Occasionally, the data was divided into groups to obtain a better 
fit. Extrapolations, whenever required, are based on the 
correlation variable. There are four modes available for costing: 

(i) Stand-alone 

(ii) Unit Operation (DOS) Reference 

(iii) User Correlation 


31 



(iv) User-supplied Cost 

For the ABE costing, the stand-alone costing option was used. 
Thus, all design parameters required by the cost models were 
specified, as were utility capacities. The calculations of 
capital investment and operating cost requires the specification 
of numerous cost parameters, or the use of ASPEN defaults for the 
same. The defaults were used in most instances. The cost indices 
used were as follows: 


Index 


Value 

Plant 


325.0 

Equipment 


346.7 

Labor 


264.9 


In addition, the labor rates assumed were ^ 13. 75/hr. 

Equipment costs for the entire ABE plant are shown in Table 
4.3. ASPEN defaults have been used to evaluate these costs, 
except for the user-specified scale-up and adjustment factors, and 
in certain cases, the material factor. Table 4.4 shows the 
installation and utility costs for the base case ABE simulation. 
As far as equipment costs are concerned, the fermentation section 
requires the most capital: about 70 percent of the total invest- 
ment is used here. The costs in this section are large mainly 
because the fermentation tanks are extremely large, and a large 
number of them are required to keep the plant in operation on a 
continuous basis. 

The major utility requirements are of cooling water and 
Steam. The largest costs (58 percent of total utility costs) are 
associated with the cooling water requirements in the fermentation 
step. Since the heat evolved in this step is of a very low 
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Table 4.3 Equipment Costs, Base Case ABE Plant 

(thousand dollars) 


EQUIPMENT 


NO. 


TOTAL PRICE 


Acid Hydrolysis: 


Crusher 

2 

341.4 

Slurry pump 

2 

105.5 

Slurry tank 

1 

43.5 

Acid hydrolyzer 

1 

72.6 

Quench tank 

1 

79.4 

Carbon column 

1 

51.8 

Centrifuge 

1 

149.6 

Pumps 

10 

65.4 

Enzyme Hydrolysis: 

Nixing tank 

1 

33.5 

Enzyme hydrolyzers 

5 

11,539.2 

Rotary filter 

1 

76.6 

Pumps 

10 

65.4 

Fermentation: 

Sterilizer 

1 

76.9 

Heat exchanger 

12 

1,090.1 

Cooling tank 

1 

. 102.4 

Permentor 

10 

69,790.3 

Pumps 

17 

355.0 

Purification: 

Column 1 

1 

191.0 

Column 2 

1 

145.4 

Column 3 

1 

299.2 

Column 4 

1 

193.4 

Column 5 

1 

299.2 

Column 6 

1 

299.2 

Evaporator 

1 

4,787.5 

Pumps 

6 

20.4 

Carbon dioxide recovery: 

Compressors 

2 

7,147.8 

Pressure vessel 

1 

261.7 

Storage: 

Solid (yard/building) 

1 

194.5 

Liquid (low pressure tanks) 

6 

3,428.0 

Gas (high pressure tanks) 

2 

- 33 - 

868.4 


909.2 


11,714.7 


71,414.7 


6,235.3 


7,409.5 


4,490.4 



Table 4.4(a) Installation Cost Summary, Base Case ABE Simulation 



(thousand 

dollars) 


SECTION 

EQUIP. COST 

MATERIAL COST 

LABOR COST 

ACID HYDROLYSIS 

909 

679 

360 

ENZYME HYDROLYSIS 

11,715 

2,065 

1,190 

FERMENTATION 

71,415 

24,092 

15,624 

PURIFICATION 

6,235 

1,393 

736 

COj RECOVERY 

7,410 

4,426 

1,960 

PRODUCT /RAW 
MATERIAL STORAGE 

4,491 

1,771 

1,071 

TESTING LABOR 

— 

- — 

2,946 

CAPITALIZED SPARES 

— 

2,403 

— 

TOTAL 

102,175 

36,829 

23,887 


Table 4.4(b) Utility Cost Summary for Base Case ABE Simulation 

(thousand dollars per year) 


SECTION 

ELECTRICITY 

STEAM 

WATER 

REPRIG 

ACID HYDROLYSIS 

207 

30,548 

644 

— 

ENZYME HYDROLYSIS 

53 

— 

2,807 

— 

FERMENTATION 

114 

11,606 

107,696 

— 

PURIFICATION 

27 

5,536 

19,439 

— 

COj RECOVERY 

— 

— 

— 

5,112 

TOTAL 

401 

47,690 

130,586 

5,112 
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quality (fermentation occurs at 33 degrees C) , the water require- 
ments are excessively large; any means of lowering these costs 
should be investigated. It should be mentioned here that it is 
probably more economical to set up a cooling tower in an actual 
plant, so that the water costs are not prohibitively large. This 
was not done in the current simulation because. of the absence of a 
module for a cooling tower in the ASPEN program. 

The capital investment report and the operating cost summary 
for this case are shown in Tables 4.5 and 4.5, which are taken 
from the actual ASPEN output report. In Table 4.5, the individual 
raw material and utility figures are not accurate, though the sum 
of the two numbers is. This is because some utilities (such as 
cooling water in the fermentors) were treated as inputs (and hence 
raw materials) to the system. Obviously, the economics do not 
change because of this. 

All costs calculated for Table 4.5 were based on ASPEN 
defaults. The hourly rate for operating labor was taken as % 
13.75, and the number of operators required per shift was esti- 
mated as 15. The depreciation method used was straight-line, with 
the plant life for tax purposes being taken as 15 years. The net 
operating cost of roughly ^253 million per year does not take 
into account the by-product credit accrued because of acetone and 
ethanol sales. If these are also taken into account, the cost of 

production of butanol works out to be 4 0.84/Kg (roughly ^ 2.55 
♦ 

per gallon), as opposed to the September 1985 selling price (20) 
of 4 0.79/Kg (about ^^.40 per gallon). There is some question as 
to whether all the by-products (such as acetic acid, butyric acid, 
acetoin, lignin, carbon-dioxide, etc.) can be sold. If it is 
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Table 4.5 Capital Investment Report^ Base Case 


M.E.T.C / VAX ASPEN SYSTEM RELEASE 6 ♦♦ PROD VERSION DATE: 12/13/83 PAGE 81 

ABE COSTING 

ECONOMIC EVALUATION SECTION 

CAPITAL INVESTMENT REPORT 
ASPEN METHOD 


• PROJECT CONSTRUCTION COSTS FOR 


LOCAT ION. I 985 QUARTER 




MATERIAL COST 


LABOR COST 

LABOR HOUR 

PROCESS UNITS 

% 

138.643.775 

% 

23.886.941 

1.737.232 

UTILITY UNITS 


0 


0 

0 

RECVNG . SHIPPING Sc STOR 


0 


0 

0 

SERVICE BUILDING 


11.390.061 


3.361.163 

244.448 

SERVICE SYST Cc DISTRBUTN 


7.714.913 


1.927.067 

140. 150 

ADDITIONAL DIRECT 


0 


0 

0 

SUBTOTAL 

% 

157.74B.751 

8 

29. 175. 172 

2. 121.831 

SITE DEVELOPMENT 


1.730.330 


34.045 

2.476 

FREIGHT 


3.334.986 




TOTAL DIRECT COST 

« 

162.814.267 

8 

29.209.217 

2, 124.307 

FIELD INDIRECT 


7.886.488 


21.322.728 

1.292.287 

FIELD CONSTR COST 
FIELD CONSTR LAB HRS 

% 

170.700.736 

8 

50.331.945 

3.416.593 

• PROJECT FIXED INVESTMENT 

FOR 

LOCATION. 

1983 QUARTER 

: i ♦* 



TOTAL FIELD CONSTRUCTION COST 
PROJECT MANAGEMENT COST » 

ENGINEERING St HOME OFFICE COST 
FEES. PERMITS St INSURANCE 
ADDITIONAL DEPRECIABLE COST 
SUBTOTAL - INDIRECTS 

14. 104.994 
28.209.988 
18.806.659 
0 

8 

221.232.700 

61*121.640 


TOTAL DEPRECIABLE (EX. CONT. ) 


8 

282.354.341 

• 

PROCESS CONTINGENCY • 

PROJECT DEFINITION CONTINGENCY 
TOTAL CONTINGENCY 

56.470.868 

28.235.434 


84.706.302 

1 

TOTAL DEPRECIABLE CAPITAL 


8 

367*060.643 

! 

1 

LAND • 

ROYALTY St EXPENSES 
ADDITIONAL NON-DEPREC I ABLE 

SUBTOTAL (NON-DEPRECIABLE COSTS) 

7.770.194 

0 

0 


7.770.194 


TOTAL FIXED INVESTMENT 


8 

374.830.B37 
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Table 4 . 6 


Annual Operating Costr Base Case 


M.E.T.C / VAX ASPEN SYSTEM 

ECONOMIC EVALUATION SECTION 
ANNUAL OFERATINS COST 
a985*» CAPACITY RATIO: 1.00) 

TOTAL RAW MATERIAL 
TOTAL UTILITIES 


RELEASE 6 PROD VERSION DATE: 12/13/8S PAGE 8 
ABE COSTING 


OPERATING LABOR ( 131490 HR) 

MAINTENANCE LABOR ( 1008173 HR) 

MAINTENANCE MATERIAL 
SUB-TOTAL MAINTENANCE 

SUPERVISION 
FRINGE BENEFITS 
OPERATING SUPPLIES 
OTHER DIRECT COST 

TOTAL DIRECT COST 

PROPERTY TAX 
PROPERTY INSURANCE 
ADM INST. It GENERAL 
DEPRECIATION (STRAIGHT LINE) 
TOTAL INDIRECT COST 

GROSS OPERATING COST 
LESS: BY-PRODUCT CREDIT 

NET OPERATING COST 


13.862,388 

9.324,233 


7,770. 194 
3, 108.078 
11.383,563 
25,900,648 


235.892.050 

25.377,712 

1.807,987 

23. 186.621 

3, 134,075 
7.521,780 
182,416 
0 

298.242.306 


48, 162.483 

346,404,789 

93.044,428 

253,360,361 


summary: 

FIXED COST 
VARIABLE COST 

TOTAL GROSS OPERATING COST 

GROSS COST EXCL. DEPRECIATION 


f 83.995.363 
262,409,426 

♦ 346.404,789 

9 320.504,141 
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assumed that only acetone and ethanol can be sold, the cost of 
production of butanol is then ^ 1.21/Kg (or 4^*^^ gallon). 
It would be reasonable to assume that the cost of production would 
be somewhere between ^ 0.84/Kg and 4 making this method 
an unattractive way of manufacturing butanol. It should be noted 
here however, that some major costs (such as cooling water and 
other utility costs) could conceivably be reduced by some better 
design concepts (such as energy integration, or the building of a 
cooling tower, etc,). The cost of production then would be 
significantly lower. 

4,6 Changes to the Basic Flowsheet 

Table 4.4 shows that the major utility cost is the cooling 
water cost in the fermentation section of the plant, and may be as 
much as 42.5 percent of the annual net operating cost (see Table 
4.6). It would therefore be obvious that any changes made to the 
flowsheet should be geared towards investigating methods of 
reducing these costs. Before that is done, it would be useful to 
see why these costs are so large. 

One advantage of continuous or semi-continuous processing in 
the case of ABE fermentations is that butanol, which is an 
inhibitor of the fermentation, may be continually removed from the 
product stream. This allows for a recycle of unconverted 
reactants, and a feed sugar concentration larger than that 
originally thought possible. There are two negative effects of 
this, however. The first is a result of the separation scheme, 
which is distillation in this case. Live steam is added in the 
first column, thereby bringing temperatures up to levels where 
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separation of the products from water may occur: this results in 
further dilution of the recycle stream, and an extra cooling cost 
because the temperature of this stream (about 100 degrees C) must 
be brought down to the fermentation temperature (33 degrees C). 
The second negative effect is that the fermentor sizes have to be 
increased, because the recycle stream is less concentrated in the 
sugars than the streams coming from the hydrolysis steps. 

One relatively easy way of handling this problem would be to 
concentrate the recycle stream by evaporating off most of the 
water. However, this would be too costly in terms of equipment 
costs as well as high pressure steam requirements. Another way 
would be to investigate possibilities of having more concentrated 
feeds to the fermentation section, by altering some operating 
conditions further upstream. Two schemes were explored, and are 
explained here. 

A. Concentrated Feed to Enzyme Hydrolyzers 

Literature surveys have shown that most enzyme hydrolysis 
studies have been carried out at a maximum Initial cellulose 
concentration of about 8 weight percent. Lee and Fan (21) mention 
that at above these concentrations, the cellulose fibers form a 
network which entraps water and renders it less movable, thereby 
hindering the hydrolysis. Particle sizes in their case was of the 
order of 0.04 mm, and the substrate involved was Solka Floe. 
Since the problem seems to be a physical one, and since our 
simulation is for roughly 2 mm particle sizes, it seems that the 
problem may not be as acute in this case. Also, the substrate is 
acid-hydrolyzed wood, and the 'netting' effect would probably not 
be seen. There has been mention in the literature (22, 23) of 
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substrate inhibition, but the data is not extensive. In any 
event, it is a worthwhile exercise to gauge the effect of having 
higher feed concentrations to the enzyme hydrolysis step on the 
entire process economics, if only to answer a 'what if' question. 
If the improvement in the economics was substantial enough, a 
potential area of research would be identified. 

As opposed to the base case of 8 wt % substrate concentration 
in the inlet to enzyme hydrolysis, a simulation with roughly 26 wt 
% was done. This had far-reaching effects on the process 
economics, ranging from different utility costs (because of lower 
steam requirements for heating during hydrolysis and for 
sterilization prior to fermentation, and of lower cooling costs 
during fermentation) to lower equipment costs (because of the 
lower throughputs). An overall view of the economics can be found 
in Tables 4.7 and 4.8, which are the capital investment report and 
the annual operating cost respectively. 

B. Two-Stage Acid Hydrolysis 

Since there is some question as to whether enzyme hydrolysis 
can proceed at high inlet substrate concentrations, an alternative 
is to carry out the hydrolysis in two stages. The amorphous 
cellulose and hemicellulose fractions could be hydrolyzed as 
before, and the remaining insoluble crystalline fraction could be 
hydrolyzed by hydrolysis at other conditions. Care can be taken 
to ensure the degradation reactions do not proceed beyond 
tolerable levels. However, the overall yield of the two-step 
process will not be as good as the combination of acid and enzyme 
hydrolysis because of some unavoidable product degradation. 
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Table 4.7 Capital Investment Report/ Case A 


M.E.T.C / VAX ASPEN SYSTEM RELEASE 6 *♦ PROD VERSION DATE: 12/13/B5 PAGE 80 

ABE COSTING 

ECONOMIC EVALUATION SECTION 

CAPITAL INVESTMENT REPORT 
ASPEN METHOD 


♦ PROJECT CONSTRUCTION COSTS FOR LOCATION, 1983 QUARTER : 1 



MATERIAL COST 

LABOR COST 

LABOR hour 

PROCESS UNITS 

UTILITY UNITS 

RECVNQ , SHIPPING STOR 

SERVICE BUILDING 

SERVICE SYST «e DISTRBUTN 

ADDITIONAL DIRECT 

f 

73,734,695 

0 

0 

6,091,274 
4, 123t847 
0 

• 13, 197,361 

0 
0 

1,797,511 

1,030,573 

0 

959,808 

0 

0 

130,728 

74,951 

0 

SUBTOTAL 

f 

83,971,816 

t 16,025,446 

1.165.487 

SITE DEVELOPMENT 
FREIGHT 


925,468 

1,783,512 

18,207 

1,324 

TOTAL DIRECT COST 

• 

B6.680.79S 

% 16,043,652 

1,166,811 

FIELD INDIRECT 


4,331,786 

11,711,866 

709,810 

FIELD CONSTR COST 
FIELD CONSTR LAB HRS 

% 

91,012,581 

% 27,755,519 

1,876,621 

• PROJECT FIXED INVESTMENT 

FOR 

LOCATION, 1985 QUARTER : 1 ♦♦ 

TOTAL FIELD CONSTRUCTION COST 
PROJECT MANAGEMENT COST 
ENGINEERING & HOME OFFICE COST 
FEES. PERMITS & INSURANCE 
ADDITIONAL DEPRECIABLE COST 
SUBTOTAL - INDIRECTS 

% 

B 

7,593,936 
15, 187,872 
10, 125,248 
0 

118.768. 100 
32.907.057 

TOTAL DEPRECIABLE <EX. 

CONT. ) 


% 

151,675,157 

PROCESS CONTINGENCY 
PROJECT DEFINITION CONTINGENCY 
TOTAL CONTINGENCY 

% 

30,335,031 
15, 167,516 

45,502,547 

TOTAL DEPRECIABLE CAPITAL 


% 

197, 177,704 

LAND 

ROYALTY & EXPENSES 
ADDITIONAL NON-DEPRECIABLE 
SUBTOTAL (NON-DEPRECIABLE 

COSTS) 

% 

4, 183,367 
0 
0 

4. 183.367 

TOTAL FIXED INVESTMENT 


B 

201,361,071 
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Tables 4,9 and 4.10 present the capital investment report and 
annual operating cost for this case. Again. there were 
improvements in the economics. These were substantial when 
compared to the base case, but not as dramatic when compared to 
Case A above. An extra cost is introduced because of the steam 
requirements for the second acid hydrolysis, but this is more than 
offset by the savings in cooling water requirements in the 
fermentation stage. The desired product yields were lower in this 
case. 

Table 4.11 compares installation costs and utility costs for 
each of the three cases. The dramatic improvement in the cost of 
the fermentation section is not matched by similar improvements in 
any other section. For Case B. there is a significant reduction 
in costs when the second hydrolysis section is. compared to the 
previous enzyme hydrolysis costs. Additionally. some real 
improvements that are not reflected in the costs are that in Case 
B. an extra species (enzyme) is not being introduced into the 
system. Therefore, enzyme recovery and production costs would be 
saved. In terms of ^/Kg. the cost of production of butanol works 
out to be as follows, assuming all the acetone and ethanol that is 
produced in each case is sold: 


Base Case 
Case A 
Case B 


with by-product 
credit 

0.84 

0.35 

0.23 


without by-product 
credit 

1.21 

0.81 

0.83 


These costs obviously do not include the effects of certain 
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Table 4.9 


Capital Investment Report, Case B 


M.E.T.C / VAX ASPEN SYSTEM RELEASE 6 ** PROD VERSION •• DATE: I2/23/B5 PAGE 75 

ASE COSTING 

ECONOMIC EVALUATION SECTION 


CAPITAL INVESTMENT REPORT 
ASPEN METHOD 

• PROJECT CONSTRUCTION COSTS FOR LOCATION* 1985 QUARTER : 1 ** 

MATERIAL COST LABOR COST LABOR HOUR 

PROCESS UNITS % 

62>063* 136 

% 11*642*432 

846*722 

UTILITY UNITS 

0 

0 

O 

RECVNG * SHIPPING «c STOP 

0 

0 

0 

SERVICE BUILDING 

5* 147* 190 

1*518*916 

110*467 

SERVICE SYST DISTRBUTN 

3*486*384 

870*845 

63*334 

ADDITIONAL DIRECT 

0 

0 

0 

SUBTOTAL S 

70*696*710 

% 14*032* 193 

1,020*523 

SITE DEVELOPMENT 
FREIGHT 

782*030 

1,507*086 

15*385 

1*119 

TOTAL DIRECT COST * 

72*985*826 

t 14*047*578 

1*021*642 

FIELD INDIRECT 

3.792.844 

10,254*732 

621*499 

FIELD CONSTR COST • 

FIELD CONSTR LAB HRS 

76*778*672 

% 24*302*311 

1*643, 141 

• PROJECT FIXED INVESTMENT FOR 

LOCATION* 1985 QUARTER : 1 

TOTAL FIELD CONSTRUCTION COST 


• 

101*080,983 

PROJECT MANAGEMENT COST 
ENGINEERING 8< HOME OFFICE COST 
FEES. PERMITS << INSURANCE 
ADDITIONAL DEPRECIABLE COST 
SUBTOTAL - INDIRECTS 

% 

6*475*819 

12.951.637 

8*634,425 

0 

28.061.880 

TOTAL DEPRECIABLE (EX. CONT. ) 


% 

129. 142.863 

PROCESS CONTINGENCY 
PROJECT DEFINITION CONTINGENCY 
TOTAL CONTINGENCY 

% 

25*828*573 

22*914*266 

38*742*859 

TOTAL DEPRECIABLE CAPITAL 


% 

167*885*722 

LAND 

ROYALTY «( EXPENSES 
ADDITIONAL NON-DEPRECIABLE 

SUBTOTAL (NON-DEPRECIABLE COSTS) 

S 

3,567*415 

0 

0 

3*567,415 

TOTAL FIXED INVESTMENT 


• 

171*453*137 
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Table 4.10 Annual Operating Cost, Case B 


M.E.T.C / VAX ASPEN SYSTEM RELEASE 6 ♦♦ PROD VERSION ♦♦ DATE: 12/23/85 PAGE 76 

ABE COSTING 

ECONOMIC EVALUATION SECTION 
ANNUAL OPERATING COST 
a?B5tt CAPACITY RATIO: 1.00) 

TOTAL RAW MATERIAL % 150«264«848 

TOTAL UTILITIES 25 » 309 1 404 


OPERATING LABOR ( 131490 HR) 

MAINTENANCE LABOR < 462867 HR) 

MAINTENANCE MATERIAL 
SUB-TOTAL MAINTENANCE 


% 6f 364* 434 

4,280*898 


1,807.987 


10,645,333 


SUPERVISION 
FRINGE BENEFITS 
OPERATING SUPPLIES 
OTHER DIRECT COST 


1*634,484 

3,922,763 

182,416 

0 


TOTAL DIRECT COST 

PROPERTY TAX 
PROPERTY INSURANCE 
ADM INST. GENERAL 

DEPRECIATION (STRAIGHT LINE) 
TOTAL INDIRECT COST 

GROSS OPERATING COST 
LESS: BY-PRODUCT CREDIT 

NET OPERATING COST 


% . 194*361,593 

♦ 3*567*415 

1*426*966 
5*936*762 
11*891*384 

22,822*526 


% 217* 184*119 

125,530*452 


• 91*653*667 


summary: 

FIXED COST 
VARIABLE COST 

TOTAL GROSS OPERATING COST 
GROSS COST EXCL. DEPRECIATION 


* 41*015*509 

176* 168*610 

t 217, 184,119 

% 205,292*735 
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Table 4.11 (a) Comparison of Installation Costs, Three Cases 

(thousand dollars) 


SECTION 

BASE CASE 

CASE A 

CASE B 

ACID HYDROLYSIS 

1948 

1948 

1948 

ENZYME HYDROLYSIS* 

14,970 

11,454 

2,155 

FERMENTATION 

111,131 

41,621 

36,407 

PURIFICATION 

8,365 

8,365 

8,365 

CO 2 RECOVERY 

13,796 

13,796 

13,796 

PRODUCT /RAW 
MATERIAL STORAGE 

7,333 

7,142 

8,843 

TESTING LABOR 

2,946 

1,576 

1,333 

CAPITALIZED SPARES 

2,403 

1,051 

860 

TOTAL 

162,892 

86,953 

73,707 


Table 4.11 (b) Comparison of Total 

(thousand 

Utility Costs, 
dollars) 

Three Cases 

SECTION 

BASE CASE 

CASE A 

CASE B 

ACID HYDROLYSIS 

31,399 

31,399 

31,399 

ENZYME HYDROLYSIS* 

2,860 

173 

10,414 

FERMENTATION 

119,416 

45,462 

33,557 

PURIFICATION 

25,002 

25,002 

25,002 

COj RECOVERY 

5,112 

5,112 

5,112 

TOTAL 

183,789 

107,148 

105,484 




For Case B, 


these are costs for the second acid hydrolysis stage 
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factors such as fermentor contamination. Also, since no standby 
equipment have been accounted for, costs would be higher. 

Talcing the selling price to be 0,79/Kg, the target cost of 
production should be around 0.69/Kg, so that a return of 15 
percent is made on investments. It would seem therefore, that 
both Case A and Case B might be viable candidates for further 
study. At the present time, both cases have associated 
difficulties. For Case A, the problem is whether the enzyme 
hydrolyzers can handle higher substrate concentrations than 8 
weight percent. For Case B, the uncertainty is in knowing whether 
the acid hydrolyzers which have worked in the bench scale will 
work in the enormous scales envisaged for this simulation. 
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APPENDIX A 


FORTRAN Listings of User Files 

This appendix contains the listings of the user files created 
to enable the simulation of the A-B-E Process. There are four 
subroutines, each written according to the guidelines in the ASPEN 
manuals. Descriptions of the input parameters are included as 
comment statements at the beginning of each subroutine. The four 
subroutines do the following: 


NAME 

. DESCRIPTION 

DSROOl 

Enzyme hydrolyzer 

USR002 

Acid hydrolyzer 

USR003 

A-B-E Fermentor 

USR004 

Sterilizer 
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nnooooooooooonnoonnonoonnnnnonnnnnn 


SUBROUTINE USROOKNS IN ,NINFI «5IN1 «5IN2 «SIN3 «SIN4 » 

1 SINFI »NSOUT ,NINFO .SOUTI . S0UT2 tS0UT3 .S0UT4 . 

2 SINFO »NSUBS t IDXSUB» ITYPE »NINT t INT ,NREAL . 

3 REAL »ID5 «NPG «NBOPST.NIW « IW «NW 

4 W «NSIZE tSIZE) 


CHANGE #1 (10/02/6S): Sign corrected on outlet heat stream 

This subroutine is the first version of an enzyme hydrolyzer 
Currently (7/12/BS) it handles only one input and one output 
stream of stream class MIXEDNC. Additionally* some require* 
are as follows: 


( 1 ) 


( 2 ) 

(3) 


(4) 

(5) 


( 6 ) 


(7) 


The component lists have to be in the following order: 

conventional: water 

NONCONVEN : crystalline cel lulose* amorphous cellu* 
lose* cel lebiose* glucose and enzyme 
Additional components may be added to these lists as 
required. 


SI units have been used. 

Enzyme activity is taken is cellulase (FPU/mL) and 
cellobiase (BU/mD* and are the first two elements of thi 
user*defined component attribute CAUSRl of NC-component 


ENZYME. 

lnt<l): no. of records to be retained for results pass 

Real(l) : lower limit for integration (hr) 

Real (2) : upper limit for integration (hr) 

Real (3) : temperature of hydrolyzer (K) 

Real (4) : pressure of hydrolyzer (Pa) 

In addition to these 4 spaces for the real vector* space 
must be kept to store profiles for the results pass. 

The additional spaces must equal at least 5*(int(l) ^ 1) 
in number. 

The dimension of RETN has to be 6<i>NCC * 37 


The kinetics for this subroutine are taken from: Pei tersen Sc 
Ross* B Sc B* 21* 997 (1979). 


Subroutines called are: GEAR* RPTHDR* USRCPY* FLASH 

C « 4 * « « • « « ^ « « « * « 4M» « « « 41 * « « « * « * * * « « « 41 * • a * « * a « a a a « « a * « « 4 4^ lie 4> « * « a * a 

IMPLICIT REALMS (A*H*0*Z> 



DIMENSION SINKl) 

*SIN2(1) *SIN3(1) 

.SlNd<l) 

1 


SINFI (1) 

*S0UT1(1) *S0UT2(1) 

« S0UT3 ( 1 ) . 

2 


S0UT4 ( 1 ) 

*SINF0(1) * IDX SUB (NSUBS) * 1 TYPE (N5UBS) * 

3 


INT(NINT) 

*REAL(NREAL) * IDS (2* 13) 

*NB0PST(3.NP0) * 

4 


IN(NIUI) 

*W(NW) *SIZE(NSIZE) 

*WS(68) 

5 


A (32) 

*XOUT(lOO) *Y(4) 

* RETN (103) 

6 


IRETN(6) 




COntlON 

/USER/ RMISS 

*IHISS *NGBAL *IPASS 

* IRESTR * ICONVG* 

1 


LMSG 

*LPMSQ *KFLAG »NHSTRY 

*NRPT *NTRMNL* 

2 


ISIZE 




COMMON 

/NCOMP/ NCC 

• NNCC 



COMMON 

/IDSCC/ IDSCC(2.1> 



COMMON 

/IDSNCC/ 1DSNCC(2>1) 



COMMON 

/IDXCC/ IDXCC(l) 



COMMON 

/IDXNCC/ IDXNCC(l) 



COMMON 

/MW/ XMW(l) 




COMMON 

/RPTGLB/ IREPFL.ISUBIXO) 
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C The next two lines are for GEAR'S benefit 

C 

COMMON /GEARLC/ XDUMB (70) « IDUMB (84) 

DATA N « TEST 1 «TE5T2« EPS* NWS «NA/4. ID- 10. IDO. . IDO. 68. 32/ 

DATA K ODE . NPK ODE . K PHASE . M X I T . ETOL . K RSLT /2 . 2 . 0 . 30 . I D-4 . 2 / 

C 

C Go to 500 for Results pass 

C 

IF (IPASS.EQ.4) GO TO 500 
L • O 

NOUT « - INT(l) 

TEMP - REAL (3) 

PRES » REAL (4) 

C 

C store initial values, and convert inlet flowrates (SI) to 

C concentrations <g/L) 

C 

HO « SINl (1)*XMW(1) 

CCO » SINl (NCC 10) 

ACO « SINl (NCC ♦ 11) 

CLO - SINl (NCC + 12) 

GCO • SINl (NCC 13) 

T * 1D3/H0 
Y(l) - CC0*T 
Y(2) » ACO^T 
Y(3) * CLO^T 
Y(4) - GCO^T 
C 

C Extract cellulase and csllobiase activities 

C 

CLLASE « SINl (NCC ^ NNCC 19) 

CLBASE e SINl (NCC f NNCC 20) 

C 

C X and XEND are the limits of the integration, passed from the 

C input file 

C 

X « REAL(l) 

XEND » REAL (2) 

JFLAG » O 
I FLAG » 0 

10 CONTINUE 

IF (IFLAQ.EQ. 1) GO TO 12 
DO 11 I - l.N 

IF (Y(I) .LE, lD-3) THEN 
TESTl - IDO 
TEST2 - ODO 
IFLAG « 1 
END IF 

11 CONTINUE 

12 CONTINUE 

CALL GEAR (N. X . XEND. Y. TESTl . TEST2.EPS. NOUT. XOUT. NWS. WS. NA. A. 

1 JFLAG. H. NSTEP. lORDER, NDERIV) 

IF (JFLAG) 900.910.20 
20 GO TO (30.40.10) JFLAG 

C 

C Evaluate the reaction rates and derivative values 

C 

30 CONTINUE 

DO 35 I - l.N 

IF (Y(I> .LE. ID-B) Yd) - ODO 
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t 


f 

t 


i 


i 


! 



33 CONTINUE 

RATEl - 0.939*Y(1)*CLLASE/(2.14 + Y(l) + 46.5*Y<2>) 

RATE2 - 3I.3^Y(2)**.6*CLLASE/(30.4 + Y(2> ♦ 9,74*Y<3)> 

RATE3 - 3.58*Y(3)»CLBASE/ (0.263 ♦ Y<3) ♦ 0.376*Y(4)) 

RATE4 » 4.12*Y(2)*CLLASE/(44. 1 + Y<2) ♦ 2.23*Y(4)) 

IF (RATEl. LE. lD-8) RATEl « ODD 

IF (RATE2.LE. lD-8) RATE2 * ODD 

IF (RATE3.LE. lD-8) RATE3 - ODO 

IF (RATE4.LE. lD-8) RATE4 « ODO 

WS(l) - - RATEl 
WS(2) « RATEl - RATE2 - RATE4 
WS(3> - RATE2 - RATE3 
WS(4) * RATE3 + RATE4 
60 TO 10 

40 TYPE •tXtY 

100 FORMAT ( IX, 6D13.3» 317) 

C 

C Stors values for results pass 

C 

REAL(5*L ♦ 5) « X 
DO 43 I - 1,4 

IF (Y(I) .LE. lD-10) Y(I) * ODO 
J * 3*L + I ♦ 3 
REAL(J) « Y(I) 

43 CONTINUE 

L « L ♦ 1 
60 TO 10 

900 WRITE (NTRMNL, 200) , JFLAG 

200 F0RMAT(1X,MNTE6RATI0N HALTED BECAUSE OF ERROR CONDITION M3) 

60 TO 999 

910 WR I TE ( NTRMNL , 300 ) 

300 FORMAT ( 1 X , » INTE6RAT ION COMPLETED SUCCESSFULLY’ ) 

C 

C Copy inlet stream into outlet stream 

C 

CALL USRCPY ( IDS ( 1 , 1 > 1 1 , 1 , NSIN, NSOUT) 

C 

C Copy the results from GEAR into the outlet stream 

C 

T « IDO/T 
S » ODO 
DO 160 I « 1,4 
Y(I) • Y(I)*T 
SOUTKNCC ♦ 1 +9) « Y(I) 

S - S Y(I) 

160 CONTINUE 

SOUTKl) ■ (CCO + ACO + CLO + 6C0 + HO - S)/XMW(1) 

N1 « NCC ♦ 1 
SOUTKNl) - ODO 
DO 170 I « 1,NCC 

SOUTKNl) * SOUTKNl) + SOUTKl) 

170 CONTINUE 

N1 « NCC ♦ 9 
N2 * N1 ♦ NNCC + 1 
S0UTKN2) - ODO 
DO 180 I > 1,NNCC 

S0UTKN2) - S0UTKN2) ♦ SOUTKNl ♦ I) 

180 CONTINUE 

SOUTKNCC + 2) * TEMP 
SOUTKNCC + 3) « PRES 
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SOUTKNCC ^ NNCC + 11) * TEMP 

SOUTKNCC ♦ NNCC + 12) « PRES 

C 

C Flash thv outlet stream to find energy requirements 

C 

CALL FLASH <S0UT1 ,NSUBS . IDXSUB. ITYPE »NBOPST.KDDE tNPKODEt 

1 KPHASE.MXIT *ETOL *TEMP ,PRES .GUESS .LMSQ . 

2 LPMSG .IRESTR.KRSLT .RETN . IRETN .LCFLAG) 

IF (LCFLAG, NE.O) WRITE (NTRMNL. 950) LCFLAG 

950 FORMAT (1X.» FLASH DID NOT CONVERGE. LCFLAG « ’.ID 

N1 * NCC ♦ NNCC 9 

SINFO(l) « -“((SOUTKNCC -► 1) ♦SOUTKNCC + 9) ♦SOUTKNCC -► 4) ♦ 

1 SOUTKNl D^SOUTKNl ♦ 4) > - (SINKNCC > !>♦ 

2 SINKNCC ♦ 9) ♦SINKNCC + 4) + SINKNl -► !)♦ 

3 SINKNl -1-4})) 

GO TO 999 

C 

C The following is the report writing section 

C 

500 CONTINUE 

II - (INT(l) -► 1) 25 

CALL RPTHDR(II.0.3.ISU6) 

WRITE (NRPT. 920) 

WRITE (NRPT. 930) . ( ( IDSNCC ( 1 . I ) . IDSNCC(2. I) ) . 1*1.4) 

920 FORMAT (/20X. ’PROFILES (G/L) ARE AS FOLLOWS:’/) 

930 FORMAT ( 8X . ’ TI ME ’ . 9X . 2A4. 6X . 2A4 . 6X . 2A4 . 6X . 2A4/ ) 

WRITE (NRPT. 940) . (REAL ( J) . J*5. 5^ ( INT ( 1 ) -U ) -►4) 

940 FORMAT (2X . 5D14. 5) 

999 RETURN 

END 
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nnnnnnnnnonnonnnnnnnnnnnnnnnnnnnnnnnnnnnnnnn 


SUBROUTINE USR002<NSIN ,N1NF1 «SINI .SIN2 «51N3 .SIN4 « 

1 SINFI tNSOUT ,NINFO .SOUTl »S0UT2 ,S0UT3 »S0UT4 , 

2 SINFO tNSUBS . IDXSUBt ITYPE «NINT « INT .NREAL « 

3 REAL .IDS «NPO .NBOPST.NIW . IW *NW 

4 W .NSIZE .SIZE) 


CHANGE MS: (9/24/6S) Avg mol wt. computed' for outlet materiel stream 
CHANGE M4: (9/24/85) H20 consumption based on total material balance 
CHANGE M3: (9/20/85) Sign corrected on outlet heat stream 
CHANGE M2: (9/20/85) GLO inserted to keep track of Mater consumption 
CHANGE Ml: (8/09/85) Kwarteng*s kinetic parameters inserted 

This subroutine is the first version of an acid hydrolyzer 
Currently (7/11/85) it handles only one input and one output 
stream of stream class tllXEDNC. Some additional requirements 
are as follows: 


( 1 ) 


( 2 ) 

(3) 


(4) 

(5) 


The component lists have to be in the following order: 
conventional: water, sulfuric acid. hmf. furfural 
NONCONVEN : crystalline cellulose, amorphous cell- 
ulose. cellobiose. glucose, enzyme, 
hemicel lulose. xylose and lignin 
Additional components may be added to these lists as 
required. 

SI units have been used 

Int(l): no. of records to be retained for results pass 

Int(2): id of substrate material 

1 : Douglas Fir 

2 : Kraft Paper 

3 : Solka Floe 

4 l Oak Sawdust 

5 : Corn Stover 

6 : Kwarteng*s kinetic constants 

Real(l) : lower limit of integration (min) 

Real (2) : i^per limit of integration (min) 

Real (3) : temperature of hydrolyzer (K) 

Real (4) : pressure of hydrolyzer (Pa) 

In addition to these 4 spaces for the real vector, space 
must be kept to store profiles for the results pass. 

The additional spaces must equal at least 7«(int(l) 1) 

in number. 


The kinetics for this routine are taken from Bhandari et al. 
BUB. 26. 320 (1984) 

See also I. K. Kwarteng. PhD Thesis. Dartmouth College (1983) 
Subroutines called: GEAR. USRCPY. FLASH. RPTHDR 


1 

2 

3 

4 

5 

6 
7 

1 


IMPLICIT REALMS (A-H. 

DIMENSION SINl (1) 
SINFI (1) 
S0UT4 ( 1 ) 
INT(NINT) 
IW(NIW) 

A (32) 

R1EXP(6) 

R2AE(6) 

COMMON /USER/ RMISS 
LMSG 


0-Z) 

.SIN2(1) 

.SOUTl (1) 

.SINFO(l) 

.REAL (NREAL) 

.W(NW) 

.XOUT(IOO) 

«R1AE(6) 

. INDEX (1) 

. IMIS5 
.LPMSG 


.SIN3(1) .5IN4(1) 

-S0UT2(1) .S0UT3(1) 

> IDXSUB(NSUBS) . ITYPE(NSUBS) 


« NGBAL 
.KFLAG 


. IDS(2.13) 

.SIZE(NSIZE) 

.Y(4) 

.R2K0(6) 

.RETN(61> 

I PASS 
NHSTRY 1 


.NB0PST(3i 
.US (68) 
«R1K0(6) 
.R2EXP(6) 

. IRETN(6) 


NPO) 


IRESTR 

NRPT 


, ICONVG. 
rNTRMNL. 
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non non#^ u nno non 


2 


ISXZE 

COMMON /NCOMP/ NCC ,NNCC 
COMMON /IDSCC/ IDSCC(2fl) 

COMMON /IDSNCC/ IDSNCC(2.l) 

COMMON /IDXCC/ IDXCC(l) 

COMMON /IDXNCC/ IDXNCCU) 

COMMON /MW/ XMW(l) 

COMMON /RPTQLB/ IREPFLf ISUB ( 10) 

COMMON /STWORK/ IDUM(14), MFt MXt MXl, MX2» MY, MCS, MNC, 

1 IDUM1<4), MIM. MIC, MIN, IDUM2(17) 

COMMON /STWKWK/ NCPM, NCPCS, NCPNC, NTRIAL. I DUM3 ( 2 ) , TC ALC , 

1 PCALC, VCALC, QCALC, BETA, RDUM(21), RESLTS(l) 

COMMON /GEARLC/ XDUMB <70) , IDUMB (84) 

EQUIVALENCE ( I NDEX < 1 ) , RESLTS (1 ) ) 

DATA K ODE , NPK ODE , K PHASE , MX I T , ETOL , KRSLT/2, 2,0,30,1 D-4 , 2/ 

DATA N,TESTl,TEST2,EPS,NWS,NA,QC/4, lD-10, IDO, . IDO, 68, 32, 1 . 9B72D-3/ 
DATA R1KO/1.73D19,2,8D20, 1,22D19,4.4D18,2.71D19, 1.45D1S/ 

DATA R2K0/2.38D14,4.9D14,3.79D14,2.8D12,2«01D14,3.84D9/ 

DATA RlEXP/1 . 34, 1 . 78, 1 . 16, 1 . 0, 2: 74 , 1 . 137/ 

DATA R2EXP/ 1 . 02, 0. 53, 0. 69 , 1 . 8, 1 . 86, 0. 569/ 

DATA RlAE/42.9,45. 1,45.2,42-9,45.3,33.72/ 

DATA R2AE/32. 87,32. 8, 32.8, 30. 0, 32. 8, 20. 99/ 

Go to 500 for Results pass 

IF (IPAS5.EQ.4) 60 TO 500 
L • O 

NOUT • • INT(l) 

TEMP « REAL (3) 

PRES « REAL (4) 

Find ths total mass flow rats 

SUM « ODO 
DO 5 1 » 1,NCC 

SUM - SUM XMWa>«SINl (I) 

CONTINUE 

DO 6 I • 1,NNCC 

SUM « SUM +SINKNCC 9 I) 

CONTINUE 

Find Msight % of acid 

ACID « 100.*SINK2)aXMW(2)/SUM 
I • INTC2) 

6CTEMP « l./(GC*TEMP) 

Evaluats ths pr s-sxponsnt ial factors of ths reactions 

RKl ■ R1K0(I)*ACID*^R1EXP(I)*DEXP< - RIAE < I ) ♦GCTEMP) 

RK2 • R2K0(I)*ACID*aR2EXP(I)aDEXP< - R2AE ( I ) *GCTEMP) 

IF (I.LE.5) THEN 

RK3 « (7.64 - 3.68/ACID)*lD20»DEXP( - 41.*GCTEMP) 

RK4 ■ <4.60 - 1.93/ACID) ♦1D14*DEXP( - 32. ♦GCTEMP) 

ELSE 

RK3 • 6.24D13^ACID^*1.1704^DEXP( - 27. 83^GCTEMP) 

RK4 • 2.33D12^ACID^^0.6S76^DEXP( - 27. 13^GCTEMP) 

END IF 

WATERF « 1./(XMW(1)4S1N1(1)) 
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C variablM 

C 

CO « SINKNCC ♦ 10)*WATERF/162. 

GO ■ (SINKNCC ♦ 11)/1A2. + SINKNCC + 13) / ISO. ) ♦WATERF 

HMFO » SINK3)*WATERF 

HO - SINKNCC ♦ 13)*WATERF/132. 

XO « SINKNCC + 16)*WATERF/130. 

FO * SINK4)*WATERF 

GLO « SINKNCC + 13 ) ♦WATERF/180. 

C 

C Ev»n though 6 species concern ust there are only 4 Independent 

C DDEs. The other 2 species (xylose and glucose) can be found 

C from a material balance 

C 

Y(l) • CO 
Y(2) • HMFO 
Y(3) » HO 
Y(4) ■ FO 
C 

C X and XEND are the limits of the integration^ passed from the 

C input file 

C 

X « REAL(l) 

XEND « REAL (2) 

JFLAG • 0 
10 CONTINUE 

DO IS I * KN 

IF (Y(I) .LE.ODO) Y(I) « ODO 
IS CONTINUE 

CALL GEAR ( N . X . XEND . Y . TEST 1 « TEST2 « EPS « NOUT . X OUT . NWS « WS . NA * A « 

1 JFLA6 . H . NSTEP 1 1 ORDER « NDER IV) 

IF (JFLAG) 900,910,20 
20 60 TO (30,40,10) JFLAG 

30 CONTINUE 

DO 3S I « 1,N 

IF (Y(I). LE.ODO) Y(I) « ODO 
3S CONTINUE 

C 

C Evaluate values of the dependent variables (glucose & xylose) 

C 

G « GO - (Y(2) - HMFO) - (Y(l) - CO) 

XX « XO - (Y(4) - FO) - (Y(3) - HO) 

C 

C Evaluate the reaction rates and derivative values 

C 

RATEl ■ RKKY(l) 

RATE2 » RK2*G 
RATE3 - RK3*Y(3) 

RATE4 « RK4«XX 
WS(1) « - RATEl 
WS(2) • RATE2 
WS(3) ■ “ RATE3 
WS(4) « RATE4 
GO TO 10 

40 G - GO - (Y(2) - HMFO) - (Y(l) - CO) 

XX - XO - (Y(4) - FO) - (Y(3) - HO) 

WRITE (NTRMNL, 100) , X, Y, Q, XX 
C 

C store values for results pass 

C 
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RCAL(7*L ♦ 5) • X 
REAL(7*L + 6) * S 
REAL<7*L ♦ 7) « XX 
DO 43 I * l»4 
J • 7*L +1+7 
REAL(J) * Y(I) 

43 CONTINUE 

L - L ♦ 1 
GO TO 10 

900 WRITE (NTRMNL. 200) t JFLA6 

GO TO 999 

910 WRITE (NTRMNLt 300) 

C 

C Copy inl»t mtrvam into outlet rntream 

C 

CALL USRCPY (IDS ( 1 • 1 ) . 1 « 1 . NSIN« NSOUT) 

C 

C Copy the results -from GEAR into the outlet stream 

C 

WATERF * l./WATERF 

SOUTKNCC ^ 10) ■ 162.*Y(1)*WATERF 

SOUTKNCC 11) *0. 

SOUTKNCC + 13) ■ 180.*Q*WATERF 
SOUTKNCC ♦ 13) * 132. »Y (3) ♦WATERF 
SOUTKNCC -► 16) « 130.*XX^WATERF 
S0UTK3) « Y (2) ♦WATERF 
S0UTK4) * Y(4) ♦WATERF 

51 » ODO 

52 « ODO 

53 » ODO 

DO 220 I « l.NCC 

IF (l.EQ. 1) GO TO 210 

52 » S2 SOUTKD^XMWd) 

53 « S3 ♦ SOUTKI) 

210 CONTINUE 

SI • SI ^ SINKI)aXhW(I) 

220 CONTINUE 

SST • S2 
II - NCC + 9 
DO 230 I « l^NNCC 

51 « SI ♦ SINKI + ID 

52 « S2 ♦ SOUTKI ♦ ID 

230 CONTINUE 

SST * SST ♦ (SI -S2) 

SOUTKI) * (SI - S2)/XI1W(1) 

S3 - S3 SOUTKI) 

SOUTKNCC ♦ 9) * SST/S3 
N1 • NCC ♦ 1 
SOUTKNl) • ODO 
DO 160 I « ItNCC 

SOUTKNl) « SOUTKNl) + SOUTKI) 

160 CONTINUE 

SOUTKNCC ♦ 2) - TEMP 
SOUTl (NCC ♦ 3) • PRES 
Nl ■ NCC ♦ 9 
N2 « Nl ♦ NNCC ♦ 1 
S0UTKN2) - ODO 
DO 170 I • l.NNCC 

S0UTKN2) « S0UTKN2) + SOUTKNl + I) 

170 CONTINUE 
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c 

c 

c 


c 

c 

c 

500 


100 

200 

300 

920 

930 

940 

950 

999 


S0UTKN2 ♦ 1) * TEMP 
S0UTKN2 2) » PRES 

Flash ths out 1st strsam to ^Ind out ensrgy rsquirsmsnts 

CALL FLASH<S0UT1 «NSUBS • XDXSUB* ITYPE «NB0PST«K0DE «NPK0DE* 

1 K PHASE, MX IT ,ET0L , TEMP ,PRES , GUESS ,LMSQ , 

2 LPMSQ , IRESTRiKRSLT ,RETN , IRETN ,LCFLAQ) 

IF (LCFLA6.NE.0) WRITE (NTRMNL, 950) , LCFLAQ 

N1 « NCC NNCC 9 

SINFO(l) « - ((SOUTKNCC ♦ D^SOUTKNCC ♦ 9)»S0UTl(NCC + 4) 

1 SOUTKNl + 1)*S0UTKN1 ♦ 4) ) - (SlNl (NCC + !)♦ 

2 SINKNCC ^ 9)«SINKNCC ^ A) * SlNl (N1 -i* 1 ) « 

3 SlNKNl ^ 4))) 

GO TO 999 


Ths ^olloNing is ths rsport writing section 
CONTINUE 

II » (INT<1) ♦ 1) + 25 
CALL RPTHDR<II,0,3, ISUB) 

WRITE (NRPT, 920) 

WRITE (NRPT, 930) 

WRITE (NRPT, 940) , (REAL ( J) , J«5, 7* ( INT ( 1 ) +1 ) ♦4) 

FORMAT (1X,F6. 2,6011.3) 

FORMAT (IX, ’INTEGRATION HALTED BECAUSE OF ERROR CONDITION M3) 
FORMATdX, MNTEGRATION COMPLETED SUCCESSFULLY’) 

FORMAT (/20X, ’PROFILES (KGMOL/L) ARE AS FOLLOWS!’/) 

FORMAT (3X. ’TIME’ ,3X,’GC0SE’ ,6X, ’ XLOSE’ ,5X, ’CCLOSE’ , 7X, 

I’HMF’ ,6X, ’HCLOSE’ ,6X, ’FURF’/) 

FORMAT ( 1 X , F6. 2, 601 1 .3) 

FORMATdX, ’FLASH DID NOT CONVERGE. LCFLAG - ’,11) 

RETURN 

END 
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OF pr/ 


SUBROUTINE USR003(NSIN «NINFI i 

1 SINFI .NSOUT »NINFO , 

2 SINFO tNSUBS .IDXSUBi 

3 REAL * IDS .NPO 

4 W .NSIZE iSlZE) 


SINl *SIN2 
SOUTl .S0UT2 
ITYPE ,NINT 
NBOPSTtNIW 


«SIN3 «8IN4 
»S0UT3 »S0UT4 
tINT «NREAL 
«IW «NW 


Thi« »ubroutin® is tbs 4irst version of tbs ABE fsrmsntor. 
Currently (8/05/B5) it hsndles only one input and two output 
streams of stream class MIXEDNC. Some additional rec^uirements 
are as follows: 




for 


only 


thi 


(1) Thi 


( 2 ) 

(3) 


SI units have 
Int (l)-Int (6) : 


(4) 


<S) 


Int (7) 

Real ( I ) 

Real (2) 

Real (3) 

Real (4) 

Real (S)-(IO) 


Real (11) 
Real (12) 
Real (13) 
Real (14) 

In addition 


to 


whose yields 
codes are: 
Product 
Acetone 
Butanol 
Ethanol 
Acetic acid 
Butyric acid 
Isopropanol 
Acetoin 
C-Dioxi de 
Hydrogen 
be retained for 


have been 


Begin Changes 8/27/65 

Real (15): now used to incorporate butanol inhibition 
Real (16) - Real (20) are now left available to the user 
whatever requirements may crop up. 

End Changes ♦♦ 

♦♦ Begin Changes 8/28/65 

N changed from 11 to 2« since other nine variables depend 
on the first. 

yin (11) inserted in dimension statement. Will accommodate 
Initial values 
♦♦ End Changes ** 

component lists have to be in the following order: 
CONVENTIONAL: water t sul fur ic acid.hmf . furfuralt 

acet one 1 butanol t ethanol « acet ic acid, 
butyric acidt i-propanol « acetoin. carbon 
dioxide .hydrogen .hydrochloric acid. 

• NONCONVEN : crystalline cel lulose. amorphous cellu- 
lose, cel lebiose. glucose. enzyme. hem i- 
cel lulose. xylose. lignin. biomass. 
Additional components may be added to these lists as 
required. 

been used. 

codes for products 
supplied. Product 
Code 
1 
2 

3 

4 

5 

6 

7 

8 
9 

no. of records 

pass 

lower limit for integration (hr) 
upper 1 imit for integrat ion (hr ) 
temperature of hydrolyzer (K) 
pressure of hydrolyzer (Pa) 
values for product yields (mol/100 mol 
glucose fermented) for 6 the products 
specified in int (l)-(6) 

% carbon recovery in products 
mu-max in Nonod growth equation 
Ks in Nonod growth equation 
m. the maintenance coefficient 
these 14 spaces for the real vector, spac 


results 


to 
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OR!<S!?^Al 

OF POOR 




PA'' 
QUAmir 


must hm kspt to store profiles for the results pass. * 
The additional spaces must equal at least 12<^ ( int (7) '4*1 ) • 
in number. m 

(6) The dimension of RETN has to be 6«NCC 37 * 

(7) In the absence of any better method» conversion of xylose 

to the products is taken to follow the same yield para- * 
meters as those estimated for fermentation of qlucose. a 

See papers by Papoutsakis et alt B Sc Bt 26t 174 (1964) and a 

B Sc Bt 27t 6B1 (1985) for a discussion on the overall equation 

used to model the fermentation process. « 

Subroutines called are: GEARt RPTHDR* USRCPYt FLASH* DGELG • 


IMPLICIT REAL*8(A-HtO-Z) 

REAL EPS2 

CHARACTER LID(l2)ai2 
DIMENSION SINl(l) tSIN2(l) 

1 SINFI(l) tSOUTKl) 

2 S0UT4(1) tSINFO(l) 

3 I NT ( N I NT > t REAL ( NRE AL ) 

4 IW(NIW) tW(NN) 

5 Adit 11) tXOUT(lOO) 

6 1RETN(6) tAAdltll) 

7 COEFF(IO) tYIELD<9) 

8 YINdl) 

COMMON /USER/ RMISS t IMISS 

1 LMSQ tLPMSQ 

2 ISIZE 

COMMON /NCOMP/ NCC tNNCC 

COMMON /IDSCC/ IDSCC(2tl) 

COMMON /IDSNCC/ IDSNCC(2tl) 
COMMON /IDXCC/ IDXCCd) 

COMMON /IDXNCC/ IDXNCC(l) 

COMMON /MW/ XMWd) 

COMMON /RPTGLB/ IREP^^Lt ISUB ( 10) 


tSIN3(l) 
t S0UT2 ( I ) 
t IDXSUB(NSUBS) 
t IDS(2f 13) 
tSIZE(NSIZE) 
tYdl) 
tBdl) 
t AG (242) 


tSlN4(l) 

t S0UT3 ( I > 

tlTYPE(NSUBS) 

tNB0PST(3tNP0) 1 

tWS(187) 

tRETN(103) 

tBB(ll) 

t IVECT(llt9) 1 


NGBAL 

KFLAQ 


t IPASS 
tNHSTRY 


t IRESTR 
tNRPT 


t ICONVGt 
tNTRMNLt 


The next two lines are for GEAR*s benefit 


1 

2 

1 

2 


COMMON /GEARLC/ XDUMB (70) t IDUMB (84) 

DATA NtTESTl,TEST2tEPSltNWStNAQ/2t ID-lOt IDOt . IDO* 187,242/ 
DATA K0DE,NPK0DEtKPHASE,MXITtET0LtKRSLT/2t2t0t30t lD-4,2/ 
DATA EPS2/1E-15/ 

DATA IVECT/7a0t lt-l,8*0dt 13*0, ItSaOt It-lt lt-lt8*0t It-lt 
12*0, 1, 12*0, 1,0, 1,5*0, 1,0, 0,2, 0,0, 1,8*0/ 

DATA LID/ ’BIOMASS* , ’GLUCOSE’ , ’ACETONE’ , ’BUTANOL’ , ’ETHANOL’ , 
’ACETIC ACID’ , ’BUTYRIC ACID’ , ’ I-PROPANOL’ , 

’ ACETOIN’ , ’C-DIOXIDE’ . ’HYDROGEN’ , ’ ATP’ / 

DATA A/1,2,ODO,2,8*ODO,-1, ItODO, 1, 10*0D0,-1,9*0D0, 1,-1, 

B*0D0,-1, 10*0D0, -1,-2, 9*0D0, -1,-2, 9*0D0,-1, 11*0D0, 

-1 , 9*0D0, -1 • -2, 8*ODO, -2, 9*0D0/ 


Go to 600 for Results pass 




4 

i 



IF (1PASS.EQ.4) GO TO 600 
L ■ 0 

NOUT « - INT(7) 

TEMP - REAL (3) 

PRES « REAL (4) 
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c 

C Set up the A matrix 

C 

DO 30 J ■ Ifll 
DO 50 I * 1*11 
AA(I.J) - A(I«J) 

50 CONTINUE 

DO 100 I « 1.6 
K - I ♦ 5 
DO 100 J « 1.11 

A(K.J) * IVECT(J.INT(D) 

AA<K«J) « A(K»J) 

100 CONTINUE 

C 

C Set up the RHS*b of the set of simultaneous eqs. 

C 

B(l) « 200./(100.-REAL(11) ) • 2. 

B(2) « O. 

B(3) - O. 

B(4) * 1.75 
B(5) * O. 

DO 150 I » 1.11 
IF (I.6T.3) THEN 

B(I) * (2.^B(1) )*REAL(I - D/IOO. 

END IF 
BB(I) « B(I) 

150 CONTINUE 

C 

C Solve the set of simultaneous linear equations 

C 

CALL D6ELG(BB«AA«11.1.EPS2.IER> 

IF (lER.NE.O) TYPE 900. 1 ER 
C 

C Calculate the coefficients of the 9 products 6 ATP 

C in the overall equation 

C 

DO 250 1-1.9 
COEFF(I) * ODO 
DO 200 K « 1 . 1 1 

COEFF(I) » COEFF(I) + I VECT <K . I ) ♦BB (K ) 

200 CONTINUE 

230 CONTINUE 

COEFF(IO) ■ 2.*BB(1) - 29.7 + BB(5) 

CFGL * 2. BB(1) 

DENOH « 1./(1B0.»CFBL) 

C 

C Evaluate the yields (g/g glucose fermented) of 

C the products 

C 

DO 300 1-1.9 

YIELD(I) « XMWd -► 4)*C0EFF(I)*DEN0M 
IF (YIELD(I).LE.ODO) YIELD(I) - ODO 
300 CONTINUE 

YBMASS - 295.0B»DENOM 
YBMINV - IDO/ YBMASS 
C 

C Initialize variables for GEAR and convert inlet 

C flowrates (SI) to concentrations (g/L) 

C 

T « 1D3/(SIN1 (1)*XMW(1)) 
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QLUC * SINKNCC 13) 

XYL * SINKNCC + 16) 

Yd) « SINKNCC + 18) ♦T 
Y<2) * <GLUC ^ XYL)*T 
DO 325 I » Kll 
IF (1.QE.3) THEN 
12 « I + 2 

Yd) - SINKI2)^XMW(I2) 

END IF 

YINd) * Yd) 

325 CONTINUE 

C 

C X «nd XEND the limits of the integrationv passed from the 

C input file 

C 

X - REALd) 

XEND « REAL (2) 

JFLAG « 0 

350 CALL GEAR < N • X * XEND « Y , TEST 1 » TEST2 > EPS 1 , NOUT • X OUT . NWS * UlS « NAG « AG « 

1 JFLAG,H,NSTEP,IORDER.NDERIV) 

IF (JFLAG) 500 » 550. 400 
400 GO TO (425.450.350) JFLAG 

C 

C Evaluate the reaction rates and derivative values 

C 

425 CONTINUE 

DELBMS « YBMINV*(Y(l) - YINd)) 

DO 430 I « N ♦ 1.11 

Yd) « YINd) + DELBMS* Y I ELD d - 2) 

430 CONTINUE 

DO 440 I - 1.11 

IF <DABS(Yd)).LE. lD-6) Yd) * ODO 
440 CONTINUE 

C BRATE - REALd2)*Y(l)*Y(2)/(REALd3) ♦ Y(2) + REAL (15) *Y (4) ) 

BRATE-d.-Y(4)/REALd5))*REALd2)*Yd)*Y<2)/(REAL(13)+Y(2) ) 
WSd) « BRATE 

WS(2) « - BRATE*YBMINV - REALd4)*Yd) 

GO TO 350 
450 CONTINUE 

C 

C store values for results pass 

C 

DELBMS « YBMINV*(Yd) - YINd)) 

DO 460 1-3.11 

Yd) • YINd) -► DELBMS*YIELDd - 2) 

460 CONTINUE 

TYPE *.X.Y 
REAL(12*L ♦ 21) * X 
DO 475 I • 1.11 
J • 12*L +1+21 

IF (DABS(Yd)).LE. lD-6) Yd) - ODO 
REAL(J) « Yd) 

475 CONTINUE 

L - L ♦ 1 
GO TO 350 

500 TYPE 910. JFLAG 

GO TO 999 
550 TYPE 920 

C 

C Copy inlet stream into outlet streams 
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1 

c 



c 

CALL USRCPYdDSdtl), 1 » 1 rNSIN, NSOUT) 

CALL USRCPY(IDS<1» 1) , 1 . 2, NSIN, NSOUT) 

C 

C Copy the results 'from GEAR into the 'first <3«seous) 

C outlet stream 

C 

T - IDO/T 

S0UTK12) - Y(10) *T/XMW(12) 

S0UT1U3) ■ Y<ll)*T/XriW(13) 

51 * ODO 

52 - ODO 

DO 560 I « IfNCC 

IF (I.LE. U.OR. I.GE. 14) SOUTKI) » ODO 

51 - SI + SOUTKI) 

52 - S2 SOUTl (I)exnW(I) 

560 CONTINUE 

SOUTKNCC + 1) ■ SI 
SOUTKNCC + 2) « TEMP 
SOUTKNCC + 3) « PRES 
SOUTKNCC P) « S2/S1 
N1 ■ NCC + 9 

DO 570 I « N1 ItNl + NNCC + 1 
SOUTKI) * ODO 
570 CONTINUE 

SOUTKNl + NNCC + 2) = TEMP 

SOUTKNl + NNCC + 3) » PRES 

C 

C Flash the 'first outlet stream 

C 

CALL FLASH (SOUTl .NSUBS . IDXSUBt ITYPE .NBOPST.KODE ,NPKODEt 

1 KPHASE^MXIT -ETOL »TEMP ,PRES .GUESS .LMSQ . 

2 LPMSQ . IRESTR.KRSLT .RETN . IRETN .LCFLAG) 

IF (LCFLAG.NE.O) TYPE 905.LCFLA6 

C 

C Copy the results from GEAR into the second outlet stream 

C 

51 « ODO 

52 « ODO 

DO 580 I - 1. 11 

IF (I.GE. 5) S0UT2(I) ■ T%Y(I - 2)/XMW(I) 

51 « SI ♦ S0UT2(I) 

52 « S2 *1* S0UT2(I)eXMW(I) 

580 CONTINUE 

DO 585 I - 13. NCC 

51 » SI -»• S0UT2(I) 

52 * S2 B0UT2(I)*XMW(I) 

585 CONTINUE 

S0UT2CNCC 1) - SI 

S0UT2(NCC ♦ 2) « TEMP 

S0UT2(NCC + 3) » PRES 

S0UT2(NCC ♦ 9) * S2/S1 

N1 » NCC -► 13 
FRAC » QLUC/(GLUC ♦ XYL) 

S0UT2(N1) ■ FRACi«Y<2)»T 

S0UT2(N1 + 3) * <1. - FRAC)*Y(2)*T 

S0UT2(N1 ♦ 5) » Y(l)eT 

S * ODO 

N1 « NCC + 9 

DO 590 I * l.NNCC 
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590 


C 

c 

c 


1 

2 


1 

2 

3 

4 


C 

C 

C 

600 


1 

2 


650 

700 


750 


S * S ♦ S0UT2(N1 -K I) 

CONTINUE 

S0UT2(N1 ♦ NNCC ♦ 1) « S 
S0UT2(N1 + NNCC -► 2) » TEMP 
S0UT2(N1 + NNCC + 3) * PRES 

Plash ths out 1st strsam and find snsrgy rsguirsmsnts 

CALL PLASH (S0UT2 »NSUBS » IDXSUBt ITYPE .NBOPSTtKODE .NPKODE. 

K PHASE, MX IT ,ETOL , TEMP ,PRES * GUESS ,LMSQ , 
LPMSG , IRESTR,KRSLT .RETN , IRETN ,LCPLAG) 

IP <LCPLAQ.NE.O) TYPE 905,LCPLA6 
Nl » NCC NNCC 9 

SINPO(l) » - (SOUTKNCC ♦ 1)*S0UT1(NCC + 9)*S0UT1(NCC -► 4) + 
S0UT2(NCC + 1)*SQUT2(NCC 9)*S0UT2<NCC ^ 4) + 
S0UT2(N1 + 1)*S0UT2(N1 ♦ 4) - (SINI (NCC 1 ) • 
SINKNCC ♦ 9)*SINl(NCC + 4) + SINI (Nl + !)• 
SINKNl -►A))) 

GO TO 999 

Ths following is ths report writing section 
CONTINUE 

J « (INT(7) + l)/5 

IP ( (5*J) .NE. (INT(7> + U) 3 ■ J ♦ I 

II - 14* J ♦ 27 

CALL RPTHDR(II,0,3,ISUB) 

CBM « 3D0 
CPGM « - CPGL 

WRITE (NRPT, 940) .LID (2) , CPGM. LID ( 1 ). CBM. YBMASS. (LID(I>. 

C0EPPCI"2).YIELD(I-2).I r 3. 1 1 ) .LID ( 12) . 
COEPP(IO) 

WRITE (NRPT. 950) 

DO 700 I « l.J 

11 « (I - 1)*60 + 21 

12 » II 48 

IP (I2.GE.NREAL) 12 « NREAL 
WR1TE(NRPT.960). (REAL(Kl) .K1 « 11.12.12) 

DO 650 K « 1,9 

11 « II 1 

12 * 12 + 1 

IP ( 12. 6E. NREAL) 12 « NREAL 

WRITE (NRPT. 970). LID (K). (REAL (K2).K2 « 11,12.12) 

CONTINUE 

CONTINUE 

11 * 21 

12 « 12*INT(7) ♦ 21 
DO 750 1-1,12 

WRITE (NTRMNL.*) . (REAL ( J) , J»1 1 , 12. 12) 

11 - II ♦ 1 

12 » 12 ♦ 1 
CONTINUE 


900 FORMAT (IX, ’WARNING IN DGELG. lER -’,13/) 

905 FORMAT (IX, ’FLASH DID NOT CONVERGE. LCPLAQ - ’,I1> 

910 PORMATdX,’ INTEGRATION HALTED BECAUSE OP ERROR CONDITION ’,13) 
920 FORMAT (IX, ’INTEGRATION COMPLETED SUCCESSFULLY’) 

930 PORMATdX, P6. 2, 2X. 11D12.3) 

940 PORMATdX, ’THE OVERALL FERMENTATION REACTION IS :’//5X, 

l’SPECIES’,20X,’C0EPP’,14X,’YIELD (G/G GLUCOSE) ’ /5X, A12, lOX, 
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2F12.3/10(SX,A12, lOX t F12. 3, 14X» D12. 3/ ) •SX.A12, 10X.F12-3/) 
930 FORMAT < IX » ’PROFILES <Q/L) ARE AS FOLLOWS:’/) 

960 F0RMAT(/1X» ’TIME’ f8X,5Dl2.4/) 

970 FORMAT aX*A12*5D12.4) 

999 RETURN 

END 
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1 

2 

3 

4 


SUBROUTINE U5R0O4(NSIN 
SINFI ,NSOUT 
SINFO tNSUBS 
REAL « IDS 
W .NSIZE 


tNINFI ,SIN1 «5IN2 
.NINFO fSOUTl fS0UT2 
•IDXSUBt ITYPE .NINT 
.NFO «NBOPST«NIW 
•SIZE) 


«SIN3 
« S0UT3 
. INT 
« IW 


«SIN4 
tS0UT4 
f NREAL 
«NW 


« 


This subroutine is the first version of an continuous steri- 
lizer. Currently (9/03/85) it handles only one inpu^ and one 
output stream of stream class MIXEDNCt and one o- infor- 

mation stream of class HEAT. Some additional requirements are 
as follows: 


( 1 ) 

( 2 ) 


(3) 

(4) 


SI units have 
ReaKl) : 

Real (2) : 

Real (3) : 

Real (4) : 

Real (5) : 

Real (6) : 

Real (7) : 

ReaKB) : 

Real (9) : 

Real (10) : 

Int(l) : 

The dimension 


been used. 

preexponential factor* microorganism death 
rate* X lD-25 i/m) 

activation energy for death rate (J/kgmol) 

temperature of sterilizer (K) 

pressure of sterilizer (Pa) 

diameter of sterilizer (m) 

diffusion coefficient (m*«2/s) 

required sterility level 

viscosity of fluid (kg/m s) 

error tolerance 

number to be subtracted from guesstimate 
of the solution of the nonlinear equation 
maximum number of Iterations for nonlinear 

search 

of RETN has to be 6«NCC 37 


Subroutines called are: USRNLE* USRCPY* FLASH 

Refer to any standard biotechnology textbook for the design 
equation for continuous sterilizers. 

^ a * a * * « « « * « * * * « e a a « a « a a s * a a « « * a a a « s a a s s s « « * « a a a a a a a ♦ > 1 . « * * * a a 41 * « 

IMPLICIT REALa8(A-H*0-Z) 

REAL EPS 

«SIN2(1) 
tSOUTKl) 
fSINFOd) 


1 

2 

3 

4 

5 

1 

2 


*SIN3(1) *9IN4(1) 

* S0UT2 ( 1 ) « S0UT3 ( 1 ) 

. IDXSUB(NSUBS) * ITYPE(NSUBS) 


* REAL (NREAL) * IDS (2* 13) 


COMMON 

COMMON 

COMMON 

COMMON 

COMWN 

COMMON 

COMMON 


.kl(NW) 

. IMISS 
. LPMSQ 

pNNCC 


*SIZE (NSIZE) 


*NB0PST(3*NP0) 
» RETN (103) 


* N6BAL 
*KFLAG 


* I PASS 
. NHSTRY 


, IRESTR 
*NRPT 


DIMENSICDN SINKl) 

SlNFl (1) 

S0UT4 ( 1 ) 

INT(NINT) 

IW(NIW) 

IRETN(6) 

COMMON /USER/ RMISS 
LMS6 
ISIZE 

/NCOMP/ NCC 
/IDSCC/ IDSCC(2*1) 

/IDSNCC/ IDSNCC(2. 1) 

/IDXCC/ IDXCC(l) 

/lOXNCC/ IDXNCCU) 

/MW/ XMW(l) 

/RPTGLB/ IREPFL.ISUB(IO) 

DATA K ODE * NPK ODE . K PHASE . MX I T , ETOL . K RSLT / 2 • 2 * O * 30 * 1 D-4 . 2/ 
DATA PI *GC/3. 1415926*8314.3/ 


ICONVG* 

NTRMNL* 


Go to 400 for results pass 

IF (IPASS.EQ.4) GO TO 500 
MAXIT ■ INT(l) 
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EPS * REAL<9) 

TEMP » REAL (3) 

PRES « REAL (4) 

DIA « REAL(S) 

C 

C Evaluate fluid valocity^ daath rata 

C 

VPLDW » SINKNCC + 1)*SIN1(NCC ♦ 9) /SINl (NCC ♦ 8) 

1 SINKNCC + NNCC ♦ 10) /SINl (NCC + NNCC * 17) 

VEL * VPL0W/(PI*DIA*^2./4.) 

Cl * REAL(2)/(GC<'TEMP) 

IF (C1.LE.50. ) THEN 

RK - REAL<1)*(DEXP(-C1)^ID25) 

ELSE 

RK « 1920-75^REAL<1)*DEXP(-C1 ♦ 50,) 

END IF 

TYPE ♦ t VEL . RK • DELTA 
C 

C Evaluate tha Raynolds numbar St diaparsion coafficiant 

C 

REYNLD - DIA*VEL«i»SINl (NCC -► B) /REAL (8) 

IF (REYNLD. LE- 2100. ) THEN 

DISP • REAL (6) + (VEL*DIA)«*2./(192.*REAL(6) ) 

ELSE 

FRICTF ■ 0.0791/ (REYNLD**. 25) 

DISP « 3. 57*VEL*DIA*DSQRT (FRICTF) 

END IF 

DELTA * DSQRTd. * 4. *RK*DISP/ (VEL**2. ) ) 

RSL « REAL (7) 

C 

C Evaluata tha initial suaas to tha solution of tha nla 

C 

XST - DL0Q(RSL*((1. * DELTA) **2. > / (4. *DELTA) ) 

XST » 2.*XST/(1. - DELTA) - REAL (10) 

C 

C Solva tha nonlinaar aquation 

C 

CALL USRNLE ( X . F , DERF t XST , EPS . MA X I T • I ER t DELTA . RSL ) 

GO TO (40,30*20) tIER + 1 
20 WRITE (NTRMNL* 900) *X*F 

TYPE 900* X*F 
GO TO 50 

30 WRITE (NTRMNL* 910) *X*F 

TYPE 910* X*F 
60 TO 50 

40 WRITE(NTRMNL*920) *X*F 

TYPE 920* X*F 
50 CONTINUE 

C 

C Calculata tha lansthm with/without disparsion 

C 

RLENQl - X*DISP/VEL 
RLENG2 - -DL0G(REAL(7))*VEL/RK 
C 

C Copy inlat straam into outlat straam 

C 

CALL USRCPY ( IDS ( 1 * 1 ) * 1 , 1 « NSIN* NSOUT) 

SOUTKNCC + 2) « TEMP 
SOUTKNCC * 3) « PRES 
SOUTKNCC ♦ NNCC * 11) « TEMP 
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c 

c 

c 


c 

c 

c 


c 

c 

c 

500 


900 

910 

920 

925 


930 


940 

950 

999 


SOUTKNCC NNCC + 12) • PRES 

Flash ths out 1st strsam to find snsrgy rsquirsmsnts 

CALL FLASHiSOUTl tNSUBS t IDXSUB* ITYPE »NB0PST»K0DE »NPK0DE» 

1 KPHASEtMXIT .ETOL •TEMP *PRES *GUaZZ .LMSG . 

2 LPMSG . IRESTR.KRSLT tRETN • IRETN iLCFLAG) 

IF (LCFLAG.NE.O) WRITE(NTRMNLt 950) LCFLAG 

N2 * NCC ^ NNCC 9 

Calculate energy requirements 

SINFO(l) * - ((SOUTKNCC ♦ 1) ♦SOUTKNCC + 9) ♦SOUTKNCC ♦ 4) + 

1 S0UTKN2 ♦ 1)*S0UTKN2 -► 4) ) - (SINKNCC ♦ !)• 

2 SINKNCC ♦ 9) ♦SINKNCC ♦ 4) ♦ SINKN2 ♦ !)♦ 

3 SINKN2 ♦ 4))> 

GO TO 999 

The following Is the report writing section 
CONTINUE 

WRITE (NRPTt 925) • (REAL(K) fK-l^B) 

WR I TE (NRPT , 930 ) • RK • VEL> REYNLD . D ISP • X • RLENQ 1 • RLENG2 
FORMAT ( IX •’ DERI V « O. XtF VALUES >>’f2D14.5) 

FORMAT (IXt ’EXCEEDED MAXIT. X»F VALUES >>’t2D14.5) 

FORMAT (IXt ’CONVERGED- X^F VALUES >>’.2D14.5) 

FORMAT ( /20X • ’KO FOR SPORE DEATH RATE «’ » 5X . D14. 5, ’ /SEC’ • 
1/20X. ’DELTA E FOR SPORE DEATH RATE •’•D14.5*’ J/KGMOL’ . /20X 
2KTEMPERATURE OF STERILIZER «’ • 3X,D14.5t ’ K ’• /20X •’ PRESSURE 
3 OF STERILIZER . 6X, D14.5. ’ PA’ • /20Xt ’ DIAMETER OF STERILIZ 

4ER «*»6X»D14.5K M’ • /20Xt ’ DIFFUSION COEFFICIENT* -’ .7XtD14. 5 
5f’ SQ M/SEC’ t/20Xt’ REQUIRED STERILITY LEVEL *’ » 4X •D14.5t /20 
6X» ’VISCOSITY OF FLUID •’ tl0XtDl4.5» ’ KQ/M SEC’ ) 

FORMAT </20X» ’SPORE DEATH RATE ■ ’ • 12X • D14- 5t ’ PER SEC’t/2 
lOXt ’FLUID VELOCITY IN STERILIZER « ’•D14.5»’ M/SEC’ • /20X • ’R 
2EYN0LDS NUMBER ■ ’• 13X* D14. 5t /20X •’ DISPERSION COEFFICIENT * 

3 ’f6X«D14.5«’ M^^2/S’ •/20Xf ’PECLET NUMBER « ’ • 15X.D14-5. /20 
4Xt ’LENGTH OF THE STERILIZER « ’.4X.D14.5*’ M WITH DISP’/20X 
5* ’LENGTH OF THE STERILIZER - ’»4X,D14.5>’ M W/OUT DISP’ > 

FORMAT (20Xt ’DISPERSION CODE «’ • 13Xt 12* /20X* ’CODE * 1 *> D 
IISPERSION PRESENT^ CODE - 2 IF ABSENT’) 

FORMAT (IXt ’FLASH DID NOT CONVERGE. LCFLAG * ’•!!) 

RETURN 

END 
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APPENDIX B 


Sample Input Data File and Output Report 

This appendix contains a sample input data file and output 
report from an ASPEN simulation. The two files for the acid 
hydrolysis simulations are presented here; similar files exist for 
the other sections. 

Since most of the components in the simulation are 
unconventional (e.g. lignin, sugars, etc.), the ASPEN data banks 
do not contain information on them. Consequently, in the input 
data file the input values for various parameters such as 
molecular weights, specific heat coefficients, etc., are defined 
in the PROP-DATA section. The property streams in the simulation 
are defined to have two substreams - one MIXED (which is of the 
same type as the streams in other process simulators), and the 
other NCPSD. The latter substream is used to represent the flow 
of nonconvent ional solids which have the particle size 
distribution attribute. The remainder of the input data file is 
similar to data files for other process simulators; if a detailed 
explanation is required, the interested reader is referred to the 
ASPEN User Manual. The output report is self-explanatory, and 
shows the material and energy balances around each unit, and the 
status of each stream in the flowsheet. 
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NEW 

TITLE ’THE ABE PROCESS’ 

DESCRIPTION ’ACID HYDROLYSIS’ 

IN-UNITS SI 
OUT-UNITS SI 

HISTORY MSG-LEVEL STREAMS ■ 6 
PROPERTIES SYS0P4 GLOBAL 

COMPONENTS CCLOSE / ACLOSE / CLBOSE / GCOSE /ENZYME /WATER H20 /& 
HCLOSE / XYLOSE / LIGNIN / H2S04 HCL / HMF CSH402 /$t 
FURFURAL C5H402 

ATTR-COMPS ENZYME CAUSRl /CCLOSE /ACLOSE /CLBOSE /GCOSE /HCLOSE /*c 
XYLOSE / LIGNIN 

COMP-NAMES CCLOSE CRYSTALLINE-CELLULOSE / ACLOSE AMORPHOUS-CELLULOSE / 
CLBOSE CELLOBIOSE / GCOSE GLUCOSE / HCLOSE HEMICELLULOSE 
NC-PROPS CCLOSE ENTHALPY ENTHLUSR / DENSITY DNSTYUSR 

NC-PROPS ACLOSE ENTHALPY ENTHLUSR / DENSITY DNSTYUSR 

NC-PROPS CLBOSE ENTHALPY ENTHLUSR / DENSITY DNSTYUSR 

NC-PROPS GCOSE ENTHALPY ENTHLUSR / DENSITY DNSTYUSR 

NC-PROPS ENZYME ENTHALPY ENTHLUSR / DENSITY DNSTYUSR 

NC-PROPS HCLOSE ENTHALPY ENTHLUSR / DENSITY DNSTYUSR 

NC-PROPS XYLOSE ENTHALPY ENTHLUSR / DENSITY DNSTYUSR 

NC-PROPS LIGNIN ENTHALPY ENTHLUSR /DENSITY DNSTYUSR 
PROP-DATA 

PROP-LIST ENTUA / DNSUA 

PVAL CCLOSE 85.0 2301.2 000/ 705.45 0000 

PVAL ACLOSE 85.0 2301.2 000/ 705.45 0000 

PVAL CLBOSE 85.0 2301.2 000/ 705.45 O O O O 

PVAL GCOSE 85.0 2301.2 000/ 705.45 0000 

PVAL HCLOSE 85.0 2301.2 000/ 705.45 0000 

PVAL XYLOSE 85.0 2301.2 000/ 705.45 O O O O 

PVAL LIGNIN 85.0 2301.2 000/ 705.45 0000 

PVAL ENZYME l.C 0.0 0 0 0 / lOOO 0000 

PROP-LIST DHFORM / DGFORM 

PVAL FURFURAL -2.092D8 / -1.255D8 
PVAL HMF -2.092D8 / -1.255D8 

PROP-LIST MW 

PVAL HMF 126.0 
PVAL H2S04 98.08 

DEF-STREAMS MIXNCPSD ALL 
DEF-SUBS-ATTR PSD PSD 

INTERVALS 10 

SI 2E-L IM ITS 0/ . 0002/ . 0004/ . 0006/ . 0008/ .001/8, 

.002/. 003/. 01/. 02/. 05 

DEF-STREAMS HEAT QHTR QONCH QACHYD QCOOLR 
DEF-STREAMS WORK SHAFT 1 


FLOWSHEET 










BLOCK 

CRUSHER 

IN 

m 

1 

OUT 

m 

2 



BLOCK 

PUMPl 

IN 

m 

2 

OUT 

m 

3 


SHAFT 1 

BLOCK 

SLRYTNK 

IN 

m 

3 4 

OUT 

m 

6 



BLOCK 

ACIDIN 

IN 

m 

6 7 

OUT 

m 

7B 



BLOCK 

HEATER 1 

IN 

m 

7B 

OUT 

m 

7C 


QHTR 

BLOCK 

ACIDHYD 

IN 

m 

7C 

OUT 

m 

e 


QACHYD 

BLOCK 

FICTSEP 

IN 

m 

e 

OUT 

m 

BA 

8B 


BLOCK 

QUENCH 

IN 

m 

BB 

OUT 

m 

9 

10 

QQNCH 

BLOCK 

COOLER 

IN 

= 

10 

OUT 

m 

lOA 

QCOOLR 

BLOCK 

CENFUQl 

IN 

as 

lOA 

OUT 

M 

11 

12 


BLOCK 

CENFUB2 

IN 

m 

12 

OUT 

m 

13 

14 


BLOCK 

CENFU63 

IN 

= 

11 13 

OUT 

m 

IS 



BLOCK 

C*-COLMN 

IN 


15 

OUT 

m 

16 

17 



STREAM 1 
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1.013D3 




SUBSTREAM MIXED TEMP * 300.13 PRES * 

MASS-FLOW WATER 45. 47 

SUBSTREAM NCPSD TEMP = 300.15 PRES » 1.013D5 

MASS-FLOW LIGNIN 7.88 / HCLOSE 10.71 /«e 

ACLOSE 7.24 / CCLOSE 19.64 

SUBS-ATTR PSD FRAC «00000000.5.5 
STREAM 4 

SUBSTREAM MIXED TEMP = 826.9415 PRES « 1 . 4D6 
MASS-FLOW WATER 38.95 

STREAM 7 

SUBSTREAM MIXED TEMP * 759.375 PRES « 1.4D6 
MASS-FLOW H2S04 1 . 9587 

BLOCK CRUSHER CRUSH 

PARAM 0.001 O 2 
BWI NCPSD 50 
BLOCK PUMPl PUMP 

PARAM PRES * 1.4D6 TYPE * 2 
BLOCK SLRYTNK MIXER 
BLOCK ACIDIN MIXER 
BLOCK HEATER 1 HEATER 

PARAM TEMP « 468.15 PRES « 1 . 4D6 
BLOCK ACIDHYD USER 

DESCRIPTION ’ACID HYDROLYZER’ 

SUBROUTINE USR002 USR002 
PARAM NINT « 2 NREAL » 46 
INT 5 6 

REAL 0.0 0.25 468.15 1.4D6 
BLOCK FICTSEP SEP 

FRAC SUBS ■ MIXED STRM « BA COMP « H2S04 FRAC * 1.0 
BLOCK QUENCH FLA5H2 

PARAM TEMP « 373.15 PRES * 1.013D5 
FRAC NCPSD 0.0 
BLOCK COOLER HEATER 

PARAM TEMP « 306.15 PRES » 1.013D5 
BLOCK CENFUQl CFUGE 

DIAMETER DIA * 1.0 RPS = 100 
CAKE-PROPS SPRES « IDIO MRES « IDIO 
RATIOS RL:R = .74 RC:R = .8 H:R * 1.0 
BLOCK CENFUG2 SEP 

FRAC SUBS«NCPSD STRM«13 COMP«CLBOSE GCOSE ENZYME XYLOSESt 
FRAC * 1.0 1.0 1.0 1.0 
BLOCK CENFUG3 MIXER 
BLOCK C-COLMN SEP 

FRAC SUBS * MIXED STRM » 16 COMP « H2S04 HMF FURFURALSt 
FRAC ■ 1.0 1.0 1.0 

STREAM-REPORT 

STANDARD OPTIONS « ALL MASS-FLOW 
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THE ABE PROCESS 
TABLE OF CONTENTS 


FLOWSHEET SECT ION 1 

FLOWSHEET CONNECTIVITY BY STREAMS 1 

FLOWSHEET CONNECTIVITY BY BLOCKS 1 

COMPUTATIONAL SEQUENCE 1 

OVERALL FLOWSHEET BALANCE 1 


PHYSICAL PROPERTIES SECTION 3 

COMPONENTS 3 

NON-CONVENT I ONAL PROPERTIES 3 

OPTION SETS 4 

UNIT OPERATIONS BLOCK SECTION 6 

CRUSHER (CRUSH > : CRUSHER 6 

PUMP (PUMP ): PUMPl 6 

MIXER (MIXER >: SLRYTNK 7 

MIXER (MIXER >; ACIDIN 6 

GENERAL-HEAT (HEATER); HEATER 1 9 

USER-MODEL (USER ) : ACIDHYD 10 

SEPARATOR (SEP ): FICTSEP 11 

flash; 2-OUTL (FLASH2); QUENCH 14 

GENERAL-HEAT (HEATER) ; COOLER IS 

CENTRIFUGE (CFUGE ): CENFUGl 16 

SEPARATOR (SEP ): CENFUG2 17 

MIXER (MIXER ): CENFUG3 20 

SEPARATOR (SEP ) : C-COLMN 21 



STREAM SECTION 2S 

DESCRIPTION OF STREAM CLASS MIXNCPSD 25 

DESCRIPTION OF STREAM CLASS HEAT 25 

DESCRIPTION OF STREAM CLASS WORK 25 

SUBSTREAM ATTR PSD TYPE: PSD 25 

2 3 6 7B 7C 26 

8 8A 8B 9 10 28 

lOA 11 12 13 14 30 

15 16 17 1 4 32 

7 34 

QHTR (PACHYD QQNCH QCOOLR 36 

SHAFTl 37 
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M.E.T.C / VAX ASPEN SYSTEM RELEASE 6 •• PROD VERSION ** DATE: 11/11/B5 PAGE 1 

THE ABE PROCESS 
FLOWSHEET SECTION 

FLOWSHEET CONNECTIVITY BY STREAMS 


STREAM 

SOURCE 

BEST 

STREAM SOURCE 

DEST 

2 

CRUSHER 

PUMPl 

3 

PUMPl 

SLRYTNK 

SHAFT 1 

PUMPl 

— 

6 

SLRYTNK 

ACIDIN 

7B 

ACIDIN 

HEATER 1 

7C 

HEATER 1 

ACIDHYD 

QHTR 

HEATER 1 

— 

8 

ACIDHYD 

FICTSEP 

QACHYD 

ACIDHYD 

— 

BA 

FICTSEP 

— — 

8B 

FICTSEP 

QUENCH 

9 

QUENCH 

— 

10 

QUENCH 

COOLER 

QQNCH 

QUENCH 

— 

lOA 

COOLER 

CENFUGl 

QCOOLR 

COOLER 

— 

11 

CENFUGl 

CENFUG3 

12 

CENFUGl 

CENFUG2 

13 

CENFUG2 

CENFUQ3 

14 

CENFUG2 

— 

15 

CENFUQ3 

C-COLMN 

16 

C-COLMN 

— 

17 

C-COLMN 

— 

1 

— 

CRUSHER 

4 

— 

SLRYTNK 

7 

— 

ACIDIN 

FLOWSHEET 

CONNECTIVITY 

BY BLOCKS 




BLOCK 

INLETS 



OUTLETS 


CRUSHER 

1 



2 


PUMPl 

2 



3 SHAFT 1 


SLRYTNK 

3 4 



6 


ACIDIN 

6 7 



7B 


HEATER 1 

7B 



7C QHTR 


ACIDHYD 

7C 



8 QACHYD 


FICTSEP 

8 



8A 8B 


QUENCH 

8B 



9 10 QQNCH 


COOLER 

10 



lOA QCOOLR 


CENFUGl 

lOA 



11 12 


CENFUQ2 

12 



13 14 


CENFUG3 

11 13 



15 


C-COLMN 

15 



16 17 


COMPUTATIONAL SEQUENCE 





SEQUENCE USED WAS: 






CRUSHER PUMPl SLRYTNK ACIDIN HEATER 1 ACIDHYD FICTSEP QUENCH COOLER 
CENFUQl CENFUQ2 CENFUQ3 C-COLMN 

OVERALL FLOWSHEET BALANCE 

•** MASS AND ENERGY BALANCE *** 


CONVENTIONAL COMPONENTS 

IN 

(KMOL/SEC) 

OUT 

RELATIVE DIFF 

WATER 

4.68609 

4.60542 

0-172165E-01 

H2S04 

0. 199704E-01 

0. 199704E-01 

O.OOOOOOE+00 

HMF 

O.OOOOOOE+00 

0.886747E-02 

-1.00000 

FURFURAL 

O.OOOOOOE+00 

0. 80025 lE-02 

-1.00000 

SUBTOTAL <K MOL/SEC ) 

4.70607 

4.64226 

0. 1355B8E-01 

(KQ/SEC ) 

86.3787 

66.8115 

-0. 498507E-02 
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M.E.T.C / VAX ASPEN SYSTEM RELEASE 6 *♦ PROD VERSION ** DATE: 11/11/03 PAGE 2 

THE ABE PROCESS 
FLOWSHEET SECTION 

OVERALL FLOWSHEET BALANCE (CONTINUED) 


NON-CONVENT I ONAL COMPONENTS (KG/SEC ) 


CCL05E 

19.6400 

17.7014 

0.987066E-01 

ACLOSE 

7.24000 

O.OOOOOOE-J-00 

1 . 00000 

CLBOSE 

0. OOOOOOE+00 

0. 000000E-»-00 

0.000000E-»-00 

GCOSE 

0. OOOOOOE-^00 

8.60230 

-1.00000 

ENZYME 

0. OOOOOOE-^OO 

0. OOOOOOE-^00 

0. OOOOOOE+00 

HCLOSE 

10.7100 

0.854619 

0. 920204 

XYLOSE 

0. OOOOOOE-^00 

9.99892 

-1.00000 

LIGNIN 

7.88000 

7.88000 

-0.2817B2E-16 

SUBTOTAL (KG/SEC ) 
TOTAL BALANCE 

43.4700 

45.0372 

0.931752E-02 

MASS (KG/SEC ) 

131.849 

131.849 

-0.293343E-09 

ENTHALPY (WATT > 

-0. 117860E+10 

-0. 117864E+10 

0.407367E-04 
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n/11/85 PAGE 6 



« 

M.E.T.C / 

VAX ASPEN 

SYSTEM RELEASE 6 PROD VERSION DATE 

THE ABE PROCESS 
UNIT OPERATIONS BLOCK SECTION 


CRUSHER 

(CRUSH 

) : CRUSHER 



INLET ■ 1 

OUTLET • 2 


PROPERTY OPTION SET SYS0P4 


MASS AND ENERGY BALANCE 


CONVENTIONAL COMPONENTS 
WATER 
H2S04 
HMF 

FURFURAL 

SUBTOTAL (KflOL/SEC) 

(KG/SEC ) 

NON-CONVENT I ONAL COMPONENTS 
CCLOSE 
ACLOSE 
CLBOSE 
GCOSE 
ENZYME 
HCLOSE 
XYLOSE 
LIGNIN 

SUBTOTAL (KG/SEC ) 

TOTAL BALANCE 

MAS5<KG/SEC ) 

ENTHALPY (WATT ) 


IN 

(KMOL/SEC) 

2.S2401 
O.OOOOOOE+00 
0. OOOOOOE+00 
O.OOOOOOE-»-00 
2. 3240 1 
45.4700 
(KG/SEC ) 
19.6400 
7.24000 
O.OOOOOOE+00 
O.OOOOOOE+00 
O.OOOOOOE+00 
10.7100 
O.OOOOOOE+00 
7. 80000 
45.4700 


90. 9400 
-0.6943B3E+09 


OUT 

2.52401 

O.OOOOOOE-i-00 

O.OOOOOOE+00 

O.OOOOOOE+OO 

2.52401 

43.4700 

19.6400 

7.24000 

O.OQOOOOE+00 

O.OOOOOOE+00 

O.OOOOOOE+OO 

10.7100 

O.OOOOOOE+OO 

7.88000 

45.4700 

90.9400 

“0.694383E+09 


RELATIVE DIFF. 

0.000000E-**00 

O.OOOOOOE-^00 

O.OOOOOOE+00 

O.OOOOOOE+00 

O.OOOOOOE+00 

0.000000E-<-00 

-0.2261 15E-16 
-0.306691E-16 
O.OOOOOOE+00 
0.000000E-»-00 
O.OOOOOOE+00 
-0.207325E-16 
O.OOOOOOE+00 
-0.2817B2E-16 
-O. 195333E-16 

O.OOOOOOE-^00 

-0.214395E-16 


♦ INPUT DATA ♦♦♦ 


OPERATING mode: 0 * PRIMARY, I = SECONDARY O 

CRUSHER TYPE: 1 * GYRATORY/ JAW, 2 « SINGLE ROLL 

3 « MULTIPLE ROLL, 4 » CAGE MILL 2 

DIAMETER OF SOLIDS OUTLET , METER 0.00100000 

BOND WORK INDEX FOR SUBSTREAM NCPSD 50.0000 


♦♦♦ RESULTS ♦♦♦ 


POWER REQUIREMENT , WATT 
PARTICLE DIAMETER WHICH IS 

LARGER THAN 807. OF INLET MASS METER 
PARTICLE DIAMETER WHICH IS 

LARGER THAN 807. OF OUTLET MASS METER 

PUMP (PUMP ) : PUMPl 

INLET • 2 OUTLET « 3 

PROPERTY OPTION SET SYS0P4 


881.439 

0.038000 

0.00086186 
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M.E.T.C / VAX ASPEN SYSTEM RELEASE 6 ♦♦ PROD VERSION DATE: 

THE ABE PROCESS 
UNIT OPERATIONS BLOCK SECTION 

PUMP (PUMP ): PUMPl (CONTINUED) 

MASS AND ENERGY BALANCE *** 



IN 

OUT 

RELATIVE DIFF 

CONVENTIONAL COMPONENTS 

(K MOL/SEC) 



WATER 

2.52401 

2.52401 

0. OOOOOOE-fOO 

H2S04 

O.OOOOOOE+00 

O.OOOOOOE-^00 

O.OOOOOOE+00 

HMF 

0. 00(X)00E+00 

O.OOOOOOE+00 

0. OOOOOOE+00 

FURFURAL 

0. OOOOOOE+00 

O.OOOOOOE+00 

0. 000000E-«-00 

SUBTOTAL ( K MOL /SEC ) 

2.52401 

2.52401 

O.OOOOOOE+00 

(KQ/SEC ) 

45.4700 

45.4700 

O.OOOOOOE+00 

NON-CONVENT I ONAL COMPONENTS (KQ/SEC ) 



CCLOSE 

19.6400 

19.6400 

O.OOOOOOE-^00 

ACLOSE 

7.24000 

7.24000 

O.OOOOOOE+00 

CLBOSE 

0. OOOOOOE-^00 

0. OOOOOOE-^00 

O.OOOOOOE+00 

GCOSE 

0. OOOOOOE^OO 

O.OOOOOOE+00 

O.OOOOOOE-^00 

ENZYME 

0. OOOOOOE-»-00 

O.OOOOOOE+00 

O.OOOOOOE+00 

HCLOSE 

10.7100 

10.7100 

O.OOOOOOE+00 

XYLOSE 

O.OOOOOOE+00 

O.OOOOOOE+00 

0. OOOOOOE+00 

LIGNIN 

7.88000 

7.88000 

O.OOOOOOE-^00 

SUBTOTAL (KQ/SEC ) 

45.4700 

45.4700 

0. OOOOOOE+00 

TOTAL BALANCE 




MASS (KQ/SEC ) 

90.9400 

90.9400 

O.OOOOOOE+00 

ENTHALPY (WATT ) 

-0. 6943B5E+09 -0. 694385E+09 

0.695742E-08 

INPUT DATA *** 



TYPE OF PUMP: 1 -CENTRIFUGAL PUMPJ 



2. SLURRY pump: 3. POSITIVE DISPLACEMENT 

PUMP 

2 

REQUIRED EXIT PRESSURE , 

N/SQM 

1,400,000. 


PUMP EFFICIENCY , 


MISSING 

DRIVER EFFICIENCY , 


1. 

ooooo 

SOLID FLOW RATE , KG/SEC 


MISSING 


RESULTS 


TYPE OF PUMP, (CAL) 

VOLUMETRIC FLOW RATE CUM/SEC 
DELTA PRESSURE N/SQM 
FLUID POWER REQUIREMENTWATT 
BRAKE POWER REQUIREMENTWATT 
ELECTRICITY REQUIREMENT WATT 
PUMP EFFICIENCY (CAL) 

MIXER (MIXER ): SLRYTNK 

INLET STREAM (S); 3 

OUTLET stream: & 

PROPERTY OPTION SET SYS0P4 


2 

O, 053706 
1,298,700. 
69,748.6 
107,305. 
107,305. 

0. 65000 
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M.E.T.C / VAX ASPEN SYSTEM RELEASE 6 PROD VERSION ♦* DATE: ll/U/BS PAGE 6 

THE ABE PROCESS 
UNIT OPERATIONS BLOCK SECTION 


MIXER (MIXER ): SLRYTNK (CONTINUED) 


♦ ♦♦ MASS AND ENERGY BALANCE 




IN 

OUT 

RELATIVE DIFF 

CONVENTIONAL COMPONENTS 

(KMOL/SEC) 



WATER 


4. 68609 

4.68609 

O.OOOOOOE+00 

H2S04 


O.OOOOOOE+00 

O.OOOOOOE+00 

O.OOOOOOE+00 

HMF 


O.OOOOOOE-^00 

O.OOOOOOE+00 

0. OOOOOOE-^00 

FURFURAL 


O.OOOOOOE+00 

O.OOOOOOE+00 

O.OOOOOOE+00 

SUBTOTAL (KMOL/SEC) 

4.68609 

4.68609 

0. OOOOOOE-^00 

(KG/SEC 

) 

84.4200 

84.4200 

O.OOOOOOE-»-00 

NON-CONVENT I ONAL 

COMPONENTS (KQ/SEC > 



CCLOSE 


19.6400 

19.6400 

O.OOOOOOE+00 

ACLOSE 


7.24000 

7.24000 

O.OOOOOOE+00 

CLBOSE 


0. OOOOOOE+00 

O.OOOOOOE+00 

O.OOOOOOE+00 

6C0SE 


0. OOOOOOE+00 

O.OOOOOOE+00 

O.OOOOOOE+00 

ENZYME 


0. OOOOOOE+00 

0. OOOOOOE+00 

O.OOOOOOE+00 

HCLOSE 


10.7100 

10.7100 

0. OOOOOOE+00 

XYLOSE 


0. OOOOOOE+00 

O.OOOOOOE+00 

O.OOOOOOE+00 

LIGNIN 


7. 88000 

7.88000 

O.OOOOOOE+00 

SUBTOTAL (KQ/SEC 

) 

43.4700 

45.4700 

0. OOOOOOE-^00 

TOTAL BALANCE 





MASS (KG/SEC ) 


129.890 

129.890 

0. OOOOOOE-^00 

ENTHALPY (WATT 

) 

-0.117691E+10 

-0. 117692E+10 

0. 197368E-05 


• INPUT DATA ♦♦♦ 

OUTLET PRESSURE ,N/SOM 
TYPE OP FLASH - TWO PHASE 
MAXIMUM NUMBER OF ITERATIONS IN FLASH 
CONVERGENCE TOLERANCE FOR FLASH 


MISSING 


30 

0.000100000 


MIXER (MIXER ): ACIDIN 

INLET STREAM (S): 6 7 

OUTLET STREAM: 7B 

PROPERTY OPTION SET SYS0P4 


MASS AND ENERGY BALANCE 


CONVENTIONAL COMPONENTS 
WATER 
H2S04 
HMF 

FURFURAL 

SUBTOTAL (KMOL/SEC) 
(KG/SEC ) 


IN 

(KMOL/SEC) 

4.68609 

0.199704E-01 

O.OOOOOOE+00 

O.OOOOOOE+00 

4.70607 

86.3787 


OUT 

4.68609 
0. 199704E-01 
O.OOOOOOE+00 
0.000000E4-00 
4.70607 
86.3787 


RELATIVE DIFF. 

0. OOOOOOE+00 
0. OOOOOOE+00 
O.OOOOOOE+00 
O.OOOOOOE+00 
O.OOOOOOE+00 
O. OOOOOOE-^00 
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MIXER <MIXER ): ACIDIN 

NON-CONVENT I ONAL COMPONENTS 
CCLOSE 
ACLOSE 
CLBOSE 
GCOSE 
ENZYME 
HCLOSE 
XYLOSE 
LIGNIN 

SUBTOTAL (KG/SEC ) 

TOTAL BALANCE 

MASS (KG/SEC ) 

ENTHALPY (WATT ) 


(CONTINUED) 
(KG/SEC > 
19.6400 
7.24000 
0. OOOOOOE+00 
0. OOOOOOE+00 
O. OOOOOOE+OO 
10.7100 
O.OOOOOOE+OO 
7.BBOOO 
43.4700 

131.849 
-0. 117849E-^10 


19.6400 
7.24000 
O.OOOOOOE+00 
O. OOOOOOE^OO 
O.OOOOOOE^OO 
10.7100 
O. OOOOOOE-^00 
7. B8000 
43.4700 

131.849 
-0. 117B49E+10 


O.OOOOOOE+00 
0. OOOOOOE-^00 
O.OOOOOOE+OO 
O.OOOOOOE+00 
O. OOOOOOE+OO 
O. OOOOOOE-**00 
oi OOOOOOE-^00 
0. OOOOOOE+OO 
0. OOOOOOE-^00 

0. OOOOOOE-^OO 
-O. 132360E-07 


INPUT DATA 

OUTLET PRESSURE .N/SQM 
TYPE OP FLASH - TWO PHASE 
MAXIMUM NUMBER OF ITERATIONS IN FLASH 
CONVERGENCE TOLERANCE FOR FLASH 


MISSING 


30 

0.000100000 


GENERAL-HEAT (HEATER) : HEATER 1 

INPUT stream: 7B 

OUTPUT stream: 7C QHTR 

PROPERTY OPTION SET SYS0P4 


HASS AND ENERGY BALANCE 




IN 

OUT 

RELATIVE DIFF 


CONVENT I ONAL COMPONENTS 

(KMOL/SEC) 




WATER 

4.68609 

4.68609 

0. OOOOOOE+00 


H2S04 

0. 199704E-01 

0. 199704E-01 

0. OOOOOOE+00 


HMF 

0. OOOOOOE-^OO 

0. 000000E-»*00 

0. OOOOOOE+OO 


FURFURAL 

0- OOOOOOE+OO 

0. OOOOOOE+00 

0. OOOOOOE^OO 

( 

r 

SUBTOTAL (KMDL/SEC) 

4. 70607 

4.70607 

0. OOOOOOE+OO 

» 

(KG/SEC ) 

86.37B7 

86.3787 

O.OOOOOOE+00 


NON-CONVENT 1 ONAL COMPONENTS (KQ/SEC ) 




CCLOSE 

19.6400 

19.6400 

0. OOOOOOE+00 

t 

ACLOSE 

7. 24000 

7.24000 

0. OOOOOOE-^OO 


CLBOSE 

0. OOOOOOE-t-00 

0. OOOOOOE4-00 

0. OOOOOOE+00 


GCOSE 

0. OOOOOOE^^OO 

0. OOOOOOE+OO 

0. OOOOOOE-^OO 


ENZYME 

O.OOOOOOE+OO 

0. OOOOOOE-i-OO 

0. OOOOOOE+00 

• 

HCLOSE 

10.7100 

10.7100 

0. OOOOOOE+00 


XYLOSE 

0. OOOOOOE+00 

0. OOOOOOE+00 

0. OOOOOOE+00 


LIGNIN 

7. 88000 

7.88000 

0. 000000E-»-00 


SUBTOTAL (KG/SEC ) 

43.4700 

43.4700 

0. OOOOOOE^OO 

f - 

TOTAL BALANCE 





MASS (KG/SEC ) 

131.849 

131.849 

0. OOOOOOE^OO 

f 

ENTHALPY (WATT ) 

-0. 117B49E+10 

-0. 117849E+10 

0.232886E-16 


i 


! 
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UNIT OPERATIONS BLOCK SECTION 


GENERAL-HEAT (HEATER): HEATER 1 (CONTINUED) 


INPUT DATA 

TWO PHASE TP FLASH 
SPECIFIED TEMPERATURE K 
SPEC I F I ED PRESSURE N/SQM 
MAXIMUM ITERATION NO. 

CONVERGENCE TOLERANCE 

TP FLASH. NO INITIAL GUESSES ARE REQUIRE&- 


468. IS 
0. 14000E+07 
30 

O. lOOOOE-03 


RESULTS 

OUTPUT TEMPERATURE K 

OUTPUT PRESSURE N/SQM 

HEAT DUTY WATT 

VAPOR FRACTION 


468. 15 
0. 14000E+07 
0. 11547E+09 
1 . 0000 


V-L PHASE EQUILIBRIUM : 

K (I) 
0.99031 
165.06 


USER-MODEL (USER ): ACIDHYD 

ACID HYDROLYZER 
INPUT STREAMS 7C 
OUTPUT STREAMS 8 QACHYD 

PROPERTY OPTION SET SYS0P4 


COMP 

WATER 

H2S04 


F(I) 

0.99576 

0.42436E-02 


X(I) 

0.99997 

0.25775E-04 


Y(I) 

0.99576 

0.42436E-02 


MASS AND ENERGY BALANCE 


CONVENTIONAL COMPONENTS 
WATER 
H2S04 
HMF 

FURFURAL 

SUBTOTAL (KMOL/SEC) 

(KQ/SEC ) 

NON-CONVENT I ONAL COMPONENTS 
CCLOSE 
ACLOSE 
CLBOSE 
GCOSE 
ENZYME 
HCLOSE 
XYLOSE 
LIGNIN 

SUBTOTAL (KG/SEC ) 

TOTAL BALANCE 

MASS (KG/SEC ) 

ENTHALPY (WATT ) 


IN 

(KMOL/SEC) 

4.68609 
0. 199704E-01 
0. OOOOOOE-*-00 
0. OOOOOOE+00 
4.70607 
86.3787 
(KG/SEC ) 
19.6400 
7.24000 
0. OOOOOOE+00 
O.OOOOOOE+00 
O.OOOOOOE-*-00 
10-7100 
O.OOOOOOE+00 
7.B8000 
45.4700 


131.849 
-0. 106302E+10 


OUT 

4.60542 
0. 199704E-01 
0.886747E-02 
0. 80025 lE-02 
4.64226 
86.8115 

17.7014 
O.OOOOOOE+00 
O.OOOOOOE-fOO 
8.60230 
0. OOOOOOE-^00 
0.854619 
9.99892 
7.88000 
45.0372 

131.849 
-0. 106302E+10 


RELATIVE DIFF. 

0. 172165E-01 
O.OOOOOOE+00 
- 1.00000 
- 1.00000 
O. 135588E-01 
-0.498507E-02 

0.987066E-01 
1 . 00000 
O.OOOOOOE^OO 
- 1.00000 
O.OOOOOOE+OO 
0.920204 
- 1.00000 
O.OOOOOOE+OO 
0.951752E-02 

O. OOOOOOE+00 
O.OOOOOOE+00 
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n.E.T.C / VAX ASPEN SYSTEM RELEASE 6 ♦♦ PROD VERSION ♦♦ DATE: 

THE ABE PROCESS 
UNIT OPERATIONS BLOCK SECTION 

U5ER--M0DEL (USER ): ACIDHYD (CONTINUED) 

PROFILES (K6M0L/L) ARE AS FOLLOWS: 


TIME 

GCOSE 

XLOSE 

CCLOSE 

HMF 

HCLOSE 

FURF 

0.00 

0.329D-03 

O-OOOD+00 

0. 144D-02 

O.OOOD-^00 

0.961D-03 

0. OOOD+OO 

0.03 

0.339D-03 

0.374D-03 

0.141D-02 

0. 204D-04 

0.5B0D-03 

0.676D-03 

0. lO 

0.347D-03 

0.388D-03 

0.138D-02 

0.412D-04 

0. 330D--03 

0. 230D-04 

0- 15 

0.354D-03 

0. 703D-03 

0.133D-02 

0.622D-04 

0.211D-03 

0. 446D-04 

0.20 

0.360D-03 

0.76SD-03 

0. 132D-02 

O.033D-O4 

0. 12BD-03 

0.690D-04 

0.23 

0.366D-03 

0. 790D-03 

0. 129D-02 

0. 103D-03 

0.767D-04 

0.948D-04 


SEPARATOR (SEP ) : FICTSEP 

INPUT STREAM - 8 
OUTPUT STREAMS - BA 
SB 

PROPERTY OPTION SET SYS0P4 


RELATIVE DIFF. 

O.OOOOOOE+00 

O.OOOOOOE^^OO 

O.OOOOOOE+00 

O.OOOOOOE+OO 


MASS AND ENERGY BALANCE 


CONVENTIONAL COMPONENTS 
WATER 
H2504 
HMF 

FURFURAL 


IN 

(KMOL/SEC) 

4.60542 
0. 199704E-01 
0. 886747E-02 
0.800231E-02 


OUT 

4.60342 
O. 199704E-01 
0.886747E-02 
0. 80023 1E-02 


1 


i 
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SEPARATOR (SEP ): FICTSEP (CONTINUED) 

NON-CONVENT I ON AL COMPONENTS (KG/SEC ) 
CCLOSE 
ACLOSE 
CLBOSE 
GCOSE 
ENZYME 
HCLOSE 
XYLOSE 
LIGNIN 

TOTAL BALANCE 

MASS (KG/SEC ) 

ENTHALPY (WATT ) 


17.7014 
O.OOOOOOE+00 
0. OOOOOOE+00 
8.60230 
O.OOOOOOE+OO 
0.85461? 
9.99892 
7.88000 

131.849 
-0. 104762E+10 


17.7014 
O.OOOOOOE+00 
O.OOOOOOE+00 
8. 60230 
0. OOOOOOE+00 
0.854619 
9.99892 
7.88000 

131.849 
-O. 104763E-^10 


O. 000000E4-00 
O. OOOOOOE-^'OO 
O.OOOOOOE+00 
0. OOOOOOE+OO 
0. 000000E-»-00 
O.OOOOOOE+OO 
O.OOOOOOE+00 
O.OOOOOOE+00 

O. OOOOOOE+00 
0.242617E-05 



♦♦♦ INPUT 

data 


FLASH SPECS FOR STREAM 

BA 



TEMPERATURE (K 

) 

MISSING 


PRESSURE (N/SQM ) 


MISSING 


VAPOR FRACTION 


MISSING 


PHASE CODE 


2 


TEMP ESTIMATE (K 

) 

MISSING 


PRES ESTIMATE (N/SQM 

) 

MISSING 


MAX. NO. ITERATIONS 


25 


CONVERGENCE TOL. 


0.000100000 


FLASH SPECS FOR STREAM 

8B 



TEMPERATURE (K 

) 

MISSING 


PRESSURE (N/SQM ) 


MISSING 


VAPOR FRACTION 


MISSING 


PHASE CODE 


2 


TEMP estimate (K 

) 

MISSING 


PRES ESTIMATE (N/SQM 

) 

MISSING 


MAX. NO. ITERATIONS 


25 


CONVERGENCE TOL. 


0.000100000 


FRACTION OF FEED 




SUBSTREAM* MIXED 




STREAM- 8A CPT- H2S04 

FRACTION- 

1 . 00000 


HEAT DUTY (WATT 


RESULTS ♦♦♦ 

2-541.72 
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THE ABE PROCESS 
UNIT OPERATIONS BLOCK SECTION 

SEPARATOR (SEP ): FICTSEP (CONTINUED) 

COMPONENT * WATER 

STREAM SUBSTREAM SPL I T FRACT I ON 

SB MIXED 1.00000 

COMPONENT « H2S04 

STREAM SUBSTREAM SPLIT FRACTION 

8A MIXED 1,00000 

COMPONENT » HMF 

STREAM SUBSTREAM SPL I T FRACT I ON 

BB MIXED 1.00000 

COMPONENT * FURFURAL 

STREAM SUBSTREAM SPLIT FRACTION 

8B MIXED 1.00000 


SPLIT FRACTION 

1 . 00000 


SPLIT FRACTION 


COMPONENT ■ CCLOSE 
STREAM SUBSTREAM 

8B NCPSD 


COMPONENT - ACLDSE 
STREAM SUBSTREAM 


COMPONENT » CLBOSE 

STREAM SUBSTREAM SPL I T FRACT I ON 
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SEPARATOR (SEP 


PICTSEP (CONTINUED) 


COMPONENT * QCOSE 

STREAM SUBSTREAM 

8B NCPSD 


SPLIT FRACTION 

1.00000 


COMPONENT 

STREAM 


ENZYME 

SUBSTREAM 


SPLIT FRACTION 


COMPONENT « HCLOSE 
STREAM SUBSTREAM 

8B NCPSD 


SPLIT FRACTION 

1.00000 


COMPONENT ■ XYLOSE 

STREAM SUBSTREAM 

8B NCPSD 


SPLIT FRACTION 

1.00000 


COMPONENT « LIGNIN 
STREAM SUBSTREAM 

8B NCPSD 


SPLIT FRACTION 

1.00000 


FLASH : 2-OUTL ( FL ASH2 > : QUENCH 

INPUT STREAM (S>: 8B 

OUTPUT STREAM (S>: 9 10 QQNCH 

PROPERTY OPTION SET SYS0P4 


*** MASS AND ENERGY BALANCE *** 


CONVENTIONAL COMPONENTS 
WATER 
H2S04 
HMF 

FURFURAL 

SUBTOTAL (KMOL/SEC) 
(K6/SEC ) 


IN 

(KMOL/SEC) 

4. 60542 
0. OOOOOOE-^00 
0.886747E-02 
0.800251E-02 
4. 62229 
84.8526 


OUT 

4.60542 
O. OOOOOOE^OO 
O. B86747E-02 
0.B00251E-02 
4. 62229 
84.8528 


RELATIVE DIFF. 

0.895404E-10 
0. 000000E-*^00 
-0.24444 I E-07 
-0.24444 I E-07 
0.240189E-16 
-0.453814E-09 
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flash: 2-OUTL <FLASH2): QUENCH (CONTINUED) 

NON-CONVENT I ONAL COMPONENTS (KG/SEC ) 


CCLOSE 

ACLOSE 

CLBOSE 

6C05E 

ENZYME 

HCLOSE 

XYLOSE 

LIGNIN 

SUBTOTAL (KG/SEC 

TOTAL BALANCE 

MASS (KG/SEC ) 
ENTHALPY (WATT 


17.7014 
0. OOOOOOE+00 
0. OOOOOOE+00 
8.60230 
0. OOOOOOE+00 
0.854619 
9.99892 
7.88000 
45.0372 

129.890 

) -0. 104387E+10 


17.7014 
O. OOOOOOE+00 
O.OOOOOOE+00 
8.60230 
O.OOOOOOE+00 
0.854619 
9-99892 
7.88000 
45.0372 

129.890 
-O- 104587E+10 


O-OOOOOOE+00 
0. 000000E-*-00 
O.OOOOOOE+00 
O.OOOOOOE+00 
O.OOOOOOE-^00 
Oi OOOOOOE+OO 
O.OOOOOOE^-00 
O.OOOOOOE+00 
O. OOOOOOE+00 

-0.297768E-09 

-0.569903E-16 


♦ INPUT DATA 

TWO PHASE TP FLASH 
SPECIFIED TEMPERATURE K 
SPECIFIED PRESSURE N/SQM 
MAXIMUM ITERATION NO. 

CONVERGENCE TOLERANCE 

TP FLASH« NO INITIAL GUESSES ARE REQUIRED. 
LIQUID ENTRAINMENT 
SOLID SPLIT FRACTIONS: 

SUBSTREAM NO. « 1 


373.15 
O. 10130E^06 
30 

O. lOOOOE-03 
O. OOOOOE+OO 


MIXED SUBSTREAM. NO SOLID SPLITS. 


SUBSTREAM NO. - 


2 vapor: o.oooooe-^oo liquid: 


0. 10000E^31 


♦♦♦ RESULTS 

OUTPUT TEMPERATURE K 

OUTPUT PRESSURE N / SQM 

HEAT DUTY WATT 

VAPOR FRACTION 


373. 15 
O- 10130E+06 
-O.2O330E+O9 
0.66936E-01 


V-L PHASE EQUILIBRIUM : 


COMP 

WATER 

HMF 

FURFURAL 


F(I) 
0.99635 
0. 191B4E-02 
O. 17313E-02 


X(I) 

0.99982 

0.94794E-04 

0.85548E-04 


Yd) 

0.94799 

0.27339E-01 

0.24672E-01 


K (I) 
0.94816 
288.39 
288.39 


GENERAL-HEAT (HEATER) : COOLER 

INPUT STREAM: 10 

OUTPUT STREAM: lOA QCOOLR 

PROPERTY OPTION SET SYS0P4 
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QENERAL-HEAT < HEATER ) : COOLER (CONTINUED) 


MASS AND ENERGY BALANCE 




IN 

OUT 

RELATIVE DIFF 

CONVENTIONAL COMPONENTS 

(KMOL/SEC) 



WATER 


4.31211 

4.31211 

O.OOOOOOE+00 

H2S04 


O.OOOOOOE+00 

O.OOOOOOE+00 

0. OOOOOOE+00 

HMF 


0. 4O8037E-O3- 

0.408837E-03 

0. OOOOOOE+00 

FURFURAL 


0.3689S7E-03 

0.368957E-03 

O.OOOOOOE^OO 

SUBTOTAL (KMOL/SEC) 

4.31289 

4.31289 

0. OOOOOOE+00 

(KG/SEC 

) 

77.7697 

77.7697 

0. OOOOOOE+00 

NON-CONVENT I ONAL 

COMPONENTS ( K 6/ SEC ) 



CCLOSE 


17.7014 

17.7014 

0. 000000E-»*00 

ACLOSE 


O.OOOOOOE+00 

O.OOOOOOE4-00 

0. OOOOOOE+00 

CLBOSE 


0. OOOOOOE+00 

0. OOOOOOE+00 

O.OOOOOOE+00 

GCOSE 


8.60230 

8.60230 

O.OOOOOOE+00 

ENZYME 


O.OOOOOOE+00 

0. OOOOOOE-^00 

0. 000000E-*-00 

HCLOSE 


0.854619 

0.854619 

0. OOOOOOE+00 

XYLOSE 


9.99892 

9.99892 

0. OOOOOOE^-00 

LIGNIN 


7.88000 

7.BB000 

O.OOOOOOE+00 

SUBTOTAL (KG/SEC 

) 

45.0372 

45.0372 

0. OOOOOOE+00 

TOTAL BALANCE 





MASS (KG/SEC ) 


122.807 

122.807 

0. OOOOOOE+00 

ENTHALPY (WATT 

) 

-0. 117575E+10 

-0. 117575E-^10 

O.OOOOOOE+00 


INPUT DATA 

TWO PHASE TP FLASH 
SPECIFIED TEMPERATURE K 
SPEC I F 1 ED PRESSURE N / SQM 

MAXIMUM ITERATION NO. 

CONVERGENCE TOLERANCE 

TP FLASH* NO INITIAL GUESSES ARE REQUIRED. 


306. 15 
0.10130E^06 
30 

0. lOOOOE-03 


RESULTS *♦* 

OUTPUT TEMPERATURE K 

OUTPUT PRESSURE N/SQM 

HEAT DUTY WATT 

VAPOR FRACTION 


306. IS 
0. 10130E*«-06 
-0.30551E+08 
O.OOOOOE+00 


V-L PHASE EQUILIBRIUM : 


COMP 

F(l> 

X(l) 

Y(l) 

K (I) 

WATER 

0.99982 

0.99982 

0.33780 

0.42790E-01 

HMF 

0.94794E-04 

0.94794E-04 

0.34808 

465.05 

FURFURAL 

0.B5S4BE-04 

0.8554BE-04 

0.31412 

465.05 


CENTR I FUQE ( CFUQE ) : CENFUQ 1 

INLET « lOA OUTLET * 11 12 

PROPERTY OPTION SET SYS0P4 
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CENTRIFUGE 


( CFUGE ) : CENFUG 1 ( CONT I NUED ) 


*** MASS AND ENERGY BALANCE ««« 




IN 

OUT 

RELATIVE DIFF 

CONVENTIONAL COMPONENTS (KMOL/SEC) 



WATER 


4.31211 

4.31211 

0. 000000E-*“00 

H2S04 


0. OOOOOOE+00 

0.000000E-*-00 

0. OOOOOOE+OO 

HMF 


0. 408837E-03 

0-40BB37E-03 

0. OOOOOOE+OO 

FURFURAL 


0. 36B957E-03 

0. 368957E-03 

0. OOOOOOE+OO 

SUBTOTAL (KMOL/SEC) 

4.31289 

4.31209 

0. OOOOOOE+OO 

(KG/SEC 

) 

77.7697 

77.7697 

-0.228413E-16 

NON-CONVENT I ONAL 

. COMPONENTS (KG/SEC ) 



CCLOSE 


17.7014 

17.7014 

0. OOOOOOE+OO 

ACLOSE 


0. OOOOOOE+00 

O.OOOOOOE+00 

0. OOOOOOE+OO 

CLBOSE 


O.OOOOOOE+00 

O.OOOOOOE'f-OO 

0. OOOOOOE+OO 

GCOSE 


8.60230 

8.60230 

0. OOOOOOE+OO 

ENZYME 


0. OOOOOOE+00 

0. OOOOOOE+OO 

0. OOOOOOE+OO 

HCLOSE 


0.854619 

0.854619 

0. OOOOOOE+OO 

XYLOSE 


9.99892 

9.99892 

0. OOOOOOE+OO 

LIGNIN 


7.88000 

7.88000 

0. OOOOOOE+OO 

SUBTOTAL (KG/SEC 

) 

45.0372 

45.0372 

0. OOOOOOE+OO 

TOTAL BALANCE 





MASS (KG/SEC 

) 

122.807 

122.807 

-0. 144646E-16 

ENTHALPY (WATT 

> -0. 120630E-*-10 

-0- 120630E-H0 

0.247055E-16 



*•* INPUT DATA ♦♦♦ 



RATIO RADIUS OF 

LIQ 

TO RADIUS OF BOW . 

0- 

74000 

RATIO-RADIUS OF 

CAKE 

TO RADIUS OP BOW . 

0. 

80000 

RATIO OF HEIGHT 

TO RADIUS OF BOW . 

1. 

00000 

CAKE RESISTANCE 

» 


0. 

100000+11 


FILTER MEDIUM RESISTANCE . 

MOISTURE CONTENT . 

POROSITY OF CAKE . 

PARTICLE SPHERICITY , 

AVERAGE PARTICLE DIAMETER t METER 

SURFACE TENSION ,N/M 

AVERAGE SOLID DENSITY .KG/CUM 

DRY SOLIDS FEED MASS FLOW RATE , KG/SEC 


O- lOOOOO+U 
MISSING 
0.45000 
0.75000 
MISSING 
0.096065 
705.450 
45.0372 


• RESULTS 


CALCULATED PARTICLE DIAMETER .METER 
RESULTED MOISTURE CONTENT . 

SELECTED BOW RADIUS . 

REVOLUTIONS PER SECOND .H2 
BASKET HEIGHT .METER 


0.00001084 
O. 0025968 
1 . 00000 
100.000 
1 . 00000 


SEPARATOR (SEP ) : CENFUG2 

INPUT STREAM - 12 
OUTPUT STREAMS - 13 
14 

PROPERTY OPTION SET SYS0P4 
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UNIT OPERATIONS BLOCK SECTION 

SEPARATOR <SEP ): CENFUQ2 (CONTINUED) 

MASS AND ENERGY BALANCE 




IN 

OUT 

RELATIVE DIFF 

CONVENTIONAL COMPONENTS 

(K MOL/SEC) 



WATER 


0.648464E-02 

0.648464E-02 

0. OOOOOOE-^00 

H2S04 


O.OOOOOOE+00 

O.OOOOOOE+00 

O.OOOOOOE-^00 

HMF 


0.614617E-06 

0.614817E-06 

O.OOOOOOE+00 

FURFURAL 


0. 554846E-06 

0.554846E-06 

0. OOOOOOE+00 

NON-CONVENT I ONAL 

COMPONENTS ( K G/SEC ) 



CCLOSE 


17.7014 

17.7014 

0. OOOOOOE-^00 

ACLOSE 


0. OOOOOOE+00 

0. OOOOOOE-fOO 

0.000000E*4*00 

CLBOSE 


0. OOOOOOE+00 

0. OOOOOOE+00 

0.00000OE-»-0O 

GCOSE 


8.60230 

8.60230 

0. OOOOOOE+00 

ENZYME 


O.OOOOOOE+00 

O.OOOOOOE-^00 

O.OOOOOOE+00 

HCLOSE 


0. B54619 

0.854619 

O.OOOOOOE+00 

XYLOSE 


9.99892 

9.99892 

O.OOOOOOE-^-OO 

LIGNIN 


7.BBOOO 

7 .80000 

O.OOOOOOE+00 

TOTAL BALANCE 





MASS (KG/SEC ) 


45. 1S42 

45. 1542 

O.OOOOOOE-*-00 

ENTHALPY (WATT 

) 

0.29B713E+08 

0.298713E^0B 

-0.623556E-16 


FLASH SPECS FOR STREAM 13 
TEMPERATURE (K ) 

PRESSURE (N/SQM ) 

VAPOR FRACTION 
PHASE CODE 

TEMP ESTIMATE (K ) 

PRES ESTIMATE (N/SCIM ) 

MAX. NO. ITERATIONS 
CONVERGENCE TOL. 

FLASH SPECS FOR STREAM 14 
TEMPERATURE (K ) 

PRESSURE (N/SQM ) 

VAPOR FRACTION 
PHASE CODE 

TEMP ESTIMATE (K ) 

PRES ESTIMATE (N/SQM ) 

MAX. NO. ITERATIONS 
CONVERGENCE TOL. 


INPUT DATA ♦♦♦ 


MISSING 

MISSING 

MISSING 

2 

MISSING 

MISSING 

25 

0.000100000 


MISSING 

MISSING 

MISSING 

2 

MISSING 

MISSING 

25 

0.000100000 
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RESULTS *** 

HEAT DUTY (WATT ) -0. 186263-08 

COMPONENT ■ WATER 

STREAM SUBSTREAM SPLIT FRACTION 

14 MIXED 1.00000 

COMPONENT ■ H2S04 

STREAM SUBSTREAM SPLIT FRACTION 

COMPONENT ■ HMF 

STREAM SUBSTREAM SPLIT FRACTION 

14 MIXED 1.00000 

COMPONENT ■ FURFURAL 

STREAM SUBSTREAM SPLIT FRACTION 

14 MIXED 1.00000 

COMPONENT ■ CCLOSE 

STREAM SUBSTREAM SPLIT FRACTION 

14 NCPSD 1.00000 

COMPONENT ■ ACLOSE 

STREAM SUBSTREAM SPLIT FRACTION 
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THE ABE. PROCESS 
UNIT OPERATIONS BLOCK SECTION 


SEPARATOR 


(SEP ) : 

CENFUQ2 

(CONTINUED) 

COMPONENT 

S 

CLBOSE 



STREAM 


SUBSTREAM 

SPLIT 

FRACTION 

COMPONENT 


GCOSE 



STREAM 


SUBSTREAM 

SPLIT 

FRACTION 

13 


NCPSO 


1 . OOOOO 

COMPONENT 

«r 

ENZYME 

* 


STREAM 


SUBSTREAM 

SPLIT 

FRACTION 

COMPONENT 

m 

HCLOSE 



STREAM 


SUBSTREAM 

SPLIT 

FRACTION 

14 


NCPSD 


1 . OOOOO 

COMPONENT 

m 

XYLOSE 



STREAM 


SUBSTREAM 

SPLIT 

FRACTION 

13 


NCPSD 


1. OOOOO 

COMPONENT 

m 

LIGNIN 



STREAM 


SUBSTREAM 

SPLIT 

FRACTION . 

14 


NCPSD 


1 . OOOOO 

tlXER 


(MIXER >: 

CENFUG3 


INLET STREAM (S): 11 


13 


OUTLET stream: 15 

PROPERTY OPTION SET SYS0P4 
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THE ABE PROCESS 
UNIT OPERATIONS BLOCK SECTION 


MIXER 


(MIXER >: CENFU63 (CONTINUED) 


*** MASS AND ENERGY BALANCE *** 




IN 

OUT 

RELATIVE DIFF 

CONVENTIONAL COMPONENTS 

(KMOL/SEC) 



WATER 


4.30563 

4.30563 

O.OOOOOOE+00 

H2S04 


O.OOOOOOE+00 

O.OOOOOOE-^00 

O.OOOOOOE+00 

HMF 


0.408222E-03 

0. 40B222E-03 

O.OOOOOOE+OO 

FURFURAL 


0.368402E-03 

0.368402E-03 

O.OOOOOOE+00 

SUBTOTAL (KMOL/SEC) 

4.30640 

4.30640 

O.OOOOOOE+00 

(KG/SEC 

) 

77.6327 

77.6327 

O.OOOOOOE+OO 

NON-CONVENT I ONAL 

COMPONENTS ( K G/SEC ) 



CCLOSE 


0. OOOOOOE+00 

O.OOOOOOE+00 

0- OOOOOOE-^00 

ACLOSE 


O.OOOOOOE+00 

O.OOOOOOE+00 

0. OOOOOOE-^OO 

CLBOSE 


O.OOOOOOE+00 

0. OOOOOOE+00 

0. OOOOOOE+OO 

GCOSE 


8.60230 

8.60230 

O-OOOOOOE+OO 

ENZYME 


0. OOOOOOE-fOO 

O. OOOOOOE+00 

O.OOOOOOE^OO 

HCLOSE 


0. OOOOOOE+00 

0. OOOOOOE+00 

O.OOOOOOE+OO 

XYLOSE 


9.99892 

9.99892 

0. OOOOOOE+OO 

LIGNIN 


O.OOOOOOE+00 

O.OOOOOOE+00 

0. OOOOOOE^OO 

SUBTOTAL (KG/SEC 

) 

18.6012 

18.6012 

0. OOOOOOE-^OO 

TOTAL BALANCE 





MASS (KG/SEC ) 


96.2539 

96. 2339 

0. OOOOOOE+OO 

ENTHALPY (WATT 

) 

-0. 122307E+10 

-0. 122307E+10 

-0.341137E-15 


*** INPUT DATA ••• 

OUTLET PRESSURE fN/SQM 
TYPE OF FLASH - TWO PHASE 
MAXIMUM NUMBER OF ITERATIONS IN FLASH 
CONVERGENCE TOLERANCE FOR FLASH 


MISSING 


30 

0.000100000 


SEPARATOR (SEP >: C-COLMN 

INPUT STREAM - 15 
OUTPUT STREAMS - 16 
17 

PROPERTY OPTION SET SYS0P4 


MASS AND ENERGY BALANCE 


IN 

CONVENTIONAL COMPONENTS (KMOL/SEC) 


WATER 

H2S04 

HMF 

FURFURAL 


4.30363 
0. OOOOOOE+00 
0. 408222E-03 
0. 368402E-03 


OUT 

4.30363 
0. OOOOOOE+00 
0.406222E-03 
0. 368402E-03 


RELATIVE DIFF. 

0. OOOOOOE+00 
O.OOOOOOE+00 
O.OOOOOOE+00 
0. OOOOOOE-^00 
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M.E.T.C / VAX ASPEN SYSTEM RELEASE 6 ♦♦ PROD VERSION ♦* DATE: 11/11/B3 PAGE 22 

THE ABE PROCESS 
UNIT OPERATIONS BLOCK SECTION 


SEPARATOR (SEP 
NON-CONVENT I ONAL 
CCLOSE 
ACLOSE 
CLBOSE 
GCOSE 
ENZYME 
HCLOSE 
XYLOSE 
LIGNIN 

TOTAL BALANCE 

MASS (KG/SEC > 
ENTHALPY (WATT 


> : C-COLMN (CONTINUED) 

COMPONENTS (KG/ SEC ) 

O.OOOOOOE+00 

O.OOOOOOE+00 

O.OOOOOOE+00 

8. 60230 
0,OOOOOOE-*-00 
O. OOOOOOE^-00 

9, 99892 

O. OOOOOOE^-00 

96.2339 

) -0. 122307E+10 


O.OOOOOOE+00 
O.OOOOOOE+00 
O. OOOOOOE-^OO 
8. 60230 
O.OOOOOOE+00 
O.OOOOOOE+OO 
9.99892 
O.OOOOOOE+OO 

96.2539 
-O. 1223UE+10 


O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O. OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 

O.OOOOOOE+OO 
O. 353066E-04 




FLASH SPECS FOR STREAM 16 
TEMPERATURE (K ) 

PRESSURE (N/SQM ) 
VAPOR FRACTION 
PHASE CODE 

TEMP ESTIMATE (K ) 

PRES ESTIMATE (N/SQM ) 

MAX. NO. ITERATIONS 
CONVERGENCE TOL. 

FLASH SPECS FOR STREAM 17 
TEMPERATURE (K ) 

PRESSURE (N/SQM > 

VAPOR FRACTION 
PHASE CODE 

TEMP ESTIMATE (K ) 

PRES ESTIMATE (N/SQM ) 

MAX. NO. ITERATIONS 
CONVERGENCE TOL. 


INPUT DATA 


MISSING 

MISSING 

MISSING 

2 

MISSING 

MISSING 

25 

0.000100000 


MISSING 

MISSING 

MISSING 

2 

MISSING 

MISSING 

25 

0.000100000 


FRACTION OF FEED 
SUBSTREAM« MIXED 
STREAM- 16 


CPT« H2S04 FR ACT I ON- 
HMF 

FURFURAL 


1 . 00000 
1.00000 
1.00000 


HEAT DUTY (WATT 


♦ RESULTS 

43.183.8 
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THE ABE PROCESS 
UNIT OPERATIONS BLOCK SECTION 

SEPARATOR < SEP ) : C-COLMN ( CONT I NUED ) 

COMPONENT ■ WATER 

STREAM SUBSTREAM SPL I T FRACT I ON 

17 MIXED 1.00000 


COMPONENT * H2S04 

STREAM SUBSTREAM 


SPLIT FRACTION 


COMPONENT « HMF 

STREAM SUBSTREAM 

16 MIXED 


SPLIT FRACTION 

I . ooooo 


COMPONENT • FURFURAL 
STREAM SUBSTREAM 

16 MIXED 


SPLIT FRACTION 

1.00000 


COMPONENT • CCLOSE 
STREAM SUBSTREAM 


SPLIT FRACTION 


COMPONENT 

STREAM 


ACLOSE 

SUBSTREAM 


SPLIT FRACTION 


COMPONENT 

STREAM 


CLBOSE 

SUBSTREAM 


SPLIT FRACTION 
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THE ABE PROCESS 
UNIT OPERATIONS BLOCK SECTION 


SEPARATOR 


(SEP ): 

C-CDLMN 

(CONTINUED) 

COMPONENT 

. 

6C0SE 



STREAM 


SUBSTREAM 

SPLIT 

FRACTION 

17 


NCPSD 


1 . 00000 

COMPONENT 

SE 

ENZYME 



STREAM 


SUBSTREAM 

SPLIT 

FRACTION 

COMPONENT 

S 

HCLOSE 



STREAM 


SUBSTREAM 

SPLIT 

FRACTION 

COMPONENT 

a 

XYLOSE 



STREAM 


SUBSTREAM 

SPLIT 

FRACTION 

17 


NCPSD 


1.00000 

COMPONENT 

« 

LIGNIN 



STREAM 


SUBSTREAM 

SPLIT 

FRACTION 
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THE ABE PROCESS 
STREAM SECTION 

DESCRIPTION OF STREAM CLASS MIXNCPSD 

STREAM CLASS : MIXNCPSD 

SUBSTREAMS : MIXED NCPSD 

SUBSTRM CLASS: MIXED NCPSD 

SUBSTRM ATTR : PSD 


DESCRIPTION OF STREAM CLASS HEAT 

STREAM CLASS : HEAT 

STREAM ATTR : HEAT 


DESCRIPTION OF STREAM CLASS WORK 

STREAM CLASS : WORK 

STREAM ATTR : WORK 



SUBSTREAM ATTR PSD TYPE: PSD 


INTERVAL 

LOWER LIMIT 

UPPER LIMIT 

1 

0,0 

METER 

.20000-03 

METER 

2 

. 20000-03 

METER 

- 40000-03 

METER 

3 

. 40000-03 

METER 

- 60000-03 

METER 

4 

. 60000-03 

METER 

- BOOOO-03 

METER 

5 

. BOOOO-03 

METER 

O.OOlO 

METER 

6 

O.OOlO 

METER 

0.0020 

METER 

7 

0,0020 

METER 

0.0030 

METER 

8 

0-0030 

METER 

0.0100 

METER 

9 

0-0100 

METER 

0.0200 

METER 

10 

0-0200 

METER 

0.0500 

METER 


I 
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THE ABE PROCESS 
STREAM SECTION 


2 3 6 7B 7C 


STREAM ID 

2 

3 

6 

7B 

7C 

FROM : 

CRUSHER 

PUMPl 

SLRYTNK 

ACIDIN 

HEATER 1 

TO : 

PUMPl 

SLRYTNK 

ACIDIN 

HEATER! 

ACIDHYD 

class: 

MIXNCPSD 

MIXNCPSD 

MIXNCPSD 

MIXNCPSD 

MIXNCPSD 


substream: mixed 

STRUCTURE 

:: CONVENTIONAL 



WATER 

KG/SEC 

45.4700 

45.4700 

84.4200 

84.4200 

84.4200 

H2S04 

KQ/SEC 

0.0 

0.0 

0.0 

1 . 9587 

1.9587 

HMF 

KG/SEC 

0.0 

0.0 

0.0 

0.0 

0.0 

FURFURAL 

KG/SEC 

0.0 

0.0 

0.0 

0.0 

0.0 

TOTAL 

KG/SEC 

45.4700 

45. 4700 

84.4200 

86. 3787 

86.3787 

TEMP 

K 

300. 1500 

300.3159 

468.2418 

467.5981 

468. 1500 

PRES 

N/SQM 

. 10130-^06 

. 14000+07 

. 14000+07 

. 14000+07 

. 14000+07 

ENTHALPY 

J/KG 

15962-^OB- 

. 15960+08-. 14522+08-. 14210+08-. 12874-^08 

VFRAC 


0.0 

0.0 

0.3233 

0.3288 

1.0000 

LFRAC 


1.0000 

1.0000 

0.6766 

0.6711 

0.0 

ENTROPY 

J/KQ-K 

-9267. 7535- 

9266.2742-5803. 1894-5655. 1047-2799.5091 

DENSITY 

KG/CUM 

846.6394 

846. 6829 

21.0611 

21.8864 

7.0786 

AVG MW 

substream: ncpsd 

18.0150 

STRUCTURE 

10.0150 18.0150 

: NON CONVENTIONAL 

18.3547 

18.3547 

CCLOSE 

KQ/SEC 

19.6400 

19.6400 

19.6400 

19.6400 

19.6400 

ACLOSE 

KG/SEC 

7.2400 

7.2400 

7.2400 

7.2400 

7.2400 

CLB05E 

KG/SEC 

0.0 

0.0 

0.0 

0.0 

0.0 

GCOSE 

KG/SEC 

0.0 

0.0 

0.0 

0.0 

0.0 

ENZYME 

KG/SEC 

0.0 

0.0 

0.0 

0.0 

0.0 

HCLOSE 

KQ/SEC 

10.7100 

10.7100 

10.7100 

10.7100 

10.7100 

XYLOSE 

KG/SEC 

0.0 

0.0 

0.0 

0.0 

0.0 

LIGNIN 

KQ/SEC 

7.8800 

7.8800 

7.8800 

7.8800 

7.8800 

TOTAL 

KG/SEC 

45.4700 

45.4700 

45.4700 

45.4700 

45.4700 

TEMP 

K 

300.1500 

300.3159 

468.2418 

467.5981 

468. 1500 

PRES 

N/SQM 

. 10130+06 

. 14000+07 

. 14000+07 

. 14000+07 

. 14000+07 

ENTHALPY 

J/KQ 

. 69079+06 

.69117+06 

. 10776+07 

.10761+07 

. 10774+07 

VFRAC 


0.0 

0.0 

0.0 

0.0 

0.0 

LFRAC 


0.0 

0.0 

0.0 

0.0 

0.0 

DENSITY 

KG/CUM 

705.4500 

705. 4500 

705.4500 

705. 4500 

705. 4500 

AVG MW 

COMPONENT ATTRIBUTES: 
ENZYME CAUSRl 

1.0000 

l.OOOO 

1.0000 

1 . 0000 

1 . 0000 

ELEMl 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM2 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM3 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM4 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM5 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM6 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM7 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM0 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM9 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEMIO 

SUBSTREAM ATTRIBUTES: 
PSD 

MISSING 

MISSING 

MISSING 

MISSING 

MISSING 


96 
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THE ABE PROCESS 
STREAM SECTION 

2 3 6 7B 7C (CONTINUED) 


STREAM ID 

2 

3 

6 

78 

7C 

FRACl 

0. 3289 

0.3289 

0.3289 

0.3289 

0.3289 

FRAC2 

0. 1451 

0. 1451 

0. 1451 

0. 1451 

0. 1451 

FRAC3 

0. 1266 

0. 1266 

0. 1266 

0. 1266 

0.1266 

FRAC4 

0. 1205 

0. 1205 

O. 1205 

O. 1205 

0. 1205 

FRAC5 

0.2542 

0.2542 

0.2542 

0.2542 

0.2542 

FRAC6 

0.0244 

0.0244 

0.0244 

0.0244 

0.0244 

FRAC7 

0.0 

0.0 

0.0 

0.0 

0.0 

FRAC8 

0.0 

0.0 

0.0 

0.0 

0.0 

FRAC9 

0.0 

0.0 

0.0 

0.0 

0.0 

FRACIO 

0.0 

0.0 

0.0 

0.0 

0.0 
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THE ABE PROCESS 
STREAM SECTION 


8 8A 8B 

9 10 






STREAM ID 

8 

BA 

88 

9 

10 

FROM : 


ACIDHYD 

FICTSEP 

FICTSEP 

QUENCH 

QUENCH 

TO : 


FICTSEP 


QUENCH 


COOLER 

class: 


MIXNCPSD 

niXNCPSD 

MIXNCPSD 

MIXNCPSD 

MIXNCPSD 

substream: mixed 

structure: 

CONVENTIONAL 



MATER 

KG/SEC 

62.966S 

0.0 

82.9665 

5.2838 

77.6826 

H2S04 

KG/SEC 

1 . 9587 

1 . 9587 

0.0 

0.0 

0.0 

HMF 

KG/SEC 

1.1173 

0.0 

1.1173 

1 . 0657 

0.0515 

FURFURAL 

KQ/SEC 

0.7688 

0.0 

0.7688 

0.7334 

0.0354 

TOTAL 

KG/SEC 

86.8114 

1 . 9587 

84.8527 

7.0830 

77.7696 

TEMP 

K 

468. 1500 

468. 1500 

468. 1500 

373. 1500 

373. 1500 

PRES 

N/SOM 

.14000-»*07 . 

14000+07 

. 14000+07 

.10130+06 

. 10130+06 

ENTHALPY 

J/KG 

12627-^OB-. 

89407+06- 

. 12898+08-. 10377+08-. 15616+08 

VFRAC 


1 . 0000 

1 . 0000 

l.OOOO 

l.OOOO 

0.0 

LFRAC 


0.0 

0.0 

0.0 

0.0 

1 . 0000 

ENTROPY 

J/KQ-K 

-2775.3276 

8.5787- 

2852 . 2458-2026 . 52 1 0-8274 . 1 375 

DENSITY 

KQ/CUM 

7.2163 

35. 9698 

7.0860 

0.7548 

800.0021 

AVG MW 


18.7002 

98. 0800 

18.3573 

22. 8932 

18.0319 


substream: ncpso 

structure: 

NON CONVENTIONAL 



CCLOSE 

KG/SEC 

17.7014 

0.0 

17.7014 

0.0 

17.7014 

ACLOSE 

KG/SEC 

0.0 

0.0 

0.0 

0.0 

0.0 

CLBOSE 

KG/SEC 

0.0 

0.0 

0.0 

0.0 

0.0 

GCOSE 

KG/SEC 

8.6022 

0.0 

8.6022 

0. 0 

8.6022 

ENZYME 

KG/SEC 

0.0 

0.0 

0.0 

0.0 

0.0 

HCLOSE 

KG/SEC 

0.8546 

0.0 

0.8546 

0.0 

0.8546 

XYLOSE 

KG/SEC 

9.9989 

0.0 

9.9989 

0.0 

9.9989 

LIGNIN 

KG/SEC 

7.8800 

0.0 

7.8800 

0.0 

7.8800 

TOTAL 

KG/SEC 

45.0372 

0.0 

45.0372 

0.0 

45. 0372 

TEMP 

K 

468. 1500 

MISSING 

468. 1500 

MISSING 

373.1500 

PRES 

N/SQM 

.14000+07 . 

14000+07 

. 14000+07 

. 10130+06 

. 10130+06 

ENTHALPY 

J/KG 

. 10774+07 

MISSING 

.10774+07 

MISSING 

. 85878+06 

VFRAC 


0.0 

MISSING 

0.0 

MISSING 

0.0 

LFRAC 


0.0 

MISSING 

0.0 

MISSING 

0.0 

DENSITY 

KG/CUM 

MISSING 

MISSING 

705.4500 

MISSING 

705.4500 

AVG MW 

COMPONENT attributes: 
ENZYME CAUSRl 

1 . 0000 

1 . 0000 

1 . 0000 

MISSING 

1 . 0000 

ELEMl 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM2 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM3 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM4 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEMS 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM6 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM7 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEMS 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM9 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEMl 0 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 


SUBSTREAM ATTRIBUTES: 
PSD 
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THE ABE PROCESS 
STREAM SECTION 


8 SA SB 9 10 (CONTINUED) 


STREAM ID 

8 

8A 

8B 

9 

10 

FRACl 

0.3289 

0.3289 

0.3289 

0.3289 

0.3289 

FRAC2 

0. 1451 

0. 1451 

0.1451 

0. 1451 

0. 1451 

FRAC3 

0. 1266 

0. 1266 

0.1266 

0. 1266 

O. 1266 

FRAC4 

0. 1203 

0. 1205 

0.1205 

0. 1205 

0. 1203 

FRACS 

0.2542 

0.2542 

0.2542 

0. 2542 

0.2542 

FRAC6 

0.0244 

0.0244 

0.0244 

0.0244 

0.0244 

FRAC7 

0.0 

0.0 

0.0 

0.0 

0.0 

FRAC8 

0.0 

0.0 

0.0 

0.0 

0.0 

FRAC9 

0.0 

0.0 

0.0 

0.0 

0.0 

FRACIO 

0.0 

0.0 

0.0 

0-0 

0.0 
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THE ABE PROCESS 
STREAM SECTION 


lOA 11 12 13 14 


STREAM ID 

lOA 

11 

12 

13 

14 

FROM : 

COOLER 

CENFUGl 

CENFUGl 

CENFUQ2 

CENFUG2 

TO : 

CENFUGl 

CENFUG3 

CENFUG2 

CENFUG3 


class: 

MIXNCPSD 

MIXNCPSD 

MIXNCPSD 

MIXNCPSD 

MIXNCPSD 


substream: mixed 

STRUCTURE 

: CONVENTIONAL 



WATER 

KG/SEC 

77.6826 

77. 5658 

0. 1168 

0.0 

0.1168 

H2S04 

KG/SEC 

0.0 

0.0 

0.0 

0.0 

0.0 

HMF 

KG/SEC 

0.0515 

0.0514 

.77467-04 

0.0 

. 77467-04 

FURFURAL 

KG/SEC 

0.0354 

0. 0353 

.53310-04 

0.0 

.53310-04 

TOTAL 

KQ/SEC 

77.7696 

77.6527 

0-1169 

0.0 

0.1169 

TEMP 

K 

306. 1500 

306. 1500 

306- 1500 

MISSING 

306. 1500 

PRES 

N/SOM 

. 10130+06 

. 10130+06 

.10130+06 

. 10130+06 

. 10130+06 

ENTHALPY 

J/KG 

-. 15919+00- 

. 15919+08-. 15919+00 

MISSING-. 15919+00 

VFRAC 


0.0 

0.0 

0.0 

MISSING 

0.0 

LFRAC 


1.0000 

1 . 0000 

1 . 0000 

MISSING 

1 . 0000 

ENTROPY 

J/KQ-K 

-9170.5223- 

9170. 5223-9170. 5223 

MISSING-9170.5223 

DENSITY 

KG/CUM 

843.4797 

843.4797 

843.4797 

MISSING 

843.4797 

AVG MW 


10.0319 

18.0319 

10.0319 

MISSING 

18.0319 


substream: ncpsd 

structure: 

NON CONVENTIONAL 



CCLOSE 

KG/SEC 

17.7014 

0.0 

17.7014 

0.0 

17.7014 

ACLOSE 

KG/SEC 

0.0 

0.0 

0.0 

0.0 

0.0 

CLBOSE 

KG/SEC 

0.0 

0.0 

0.0 

0-0 

0.0 

GCDSE 

KG/SEC 

8.6022 

0.0 

8.6022 

8.6022 

0.0 

ENZYME 

KG/SEC 

0.0 

0.0 

0.0 

0.0 

0.0 

HCLOSE 

KG/SEC 

0.8546 

0.0 

0.8546 

0.0 

0.8546 

XYLOSE 

KG/SEC 

9.9909 

0.0 

9.9909 

9.9989 

0.0 

LIGNIN 

KG/SEC 

7.8000 

0.0 

7.0800 

0.0 

7.0800 

TOTAL 

KG/SEC 

45.0372 

0-0 

45.0372 

10.6012 

26.4360 

TEMP 

K 

306. 1500 

MISSING 

306. 1500 

306. 1500 

306. 1500 

PRES 

N/SQM 

. 10130+06 . 

10130+06 

. 10130+06 

. 10130+06 

. 10130+06 

ENTHALPY 

J/KQ 

. 70460+06 

MISSING 

.70460+06 

. 70460+06 

. 70460+06 

VFRAC 


0.0 

MISSING 

0.0 

0.0 

0.0 

LFRAC 


0.0 

MISSING 

0.0 

0.0 

0.0 

DENSITY 

KG/CUM 

705.4500 

MISSING 

705.4500 

705. 4500 

705.4500 

AVG MM 

COMPONENT attributes: 
ENZYME CAUSRl 

1.0000 

1.0000 

1 . OOOO 

1 . OOOO 

1 . OOOO 

ELEMl 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM2 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM3 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM4 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM5 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM6 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM7 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM8 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM9 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEMIO 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 


SUBSTREAM ATTRIBUTES: 
PSD 
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THE ABE PROCESS 
STREAM SECTION 


lOA 11 12 13 14 (CONTINUED) 


STREAM ID 

lOA 

11 

12 

13 

14 

PRACl 

0.3289 

MISSING 

0.3289 

0.3289 

0.3289 

FRAC2 

0. 1451 

MISSING 

0. 1451 

0. 1451 

0. 1431 

PRAC3 

0. 1266 

MISSING 

0. 1266 

0. 1266 

0. 1266 

FRAC4 

0. 1205 

MISSING 

0. 1205 

0. 1203 

0.1205 

FRACS 

0.2342 

MISSING 

0.2542 

0.2542 

0.2542 

FRAC6 

0.0244 

MISSING 

0.0244 

0.0244 

0.0244 

FRAC7 

0.0 

MISSING 

0.0 

0.0 

0.0 

FRACS 

0.0 

MISSING 

0.0 

0.0 

0.0 

FRAC9 

0.0 

MISSING 

0.0 

0.0 

0.0 

FRACIO 

0.0 

MISSING 

0.0 

0.0 

0.0 


101 



ri.E.T.C / VAX ASPEN SYSTEM RELEASE 6 ** PROD VERSION •• DATE! 11 / 11/83 PAGE 32 

THE ABE PROCESS 
STREAM SECTION 

IS 16 17 1 4 


STREAM ID 

15 

16 

17 

1 

4 

FROM : 

CENFUG3 

C-COLMN 

C-COLMN 



TO : 

C-COLMN 



CRUSHER 

SLRYTNK 

class: 

MIXNCPSD 

MIXNCPSD 

MIXNCPSD 

MIXNCPSD 

MIXNCPSD 


substream: mixed 

STRUCTURE 

! CONVENTIONAL 



water 

KG/SEC 

77.5658 

0.0 

77.5658 

45.4700 

38.9500 

H2S04 

KG/SEC 

0.0 

0.0 

0.0 

0.0 

0.0 

HMF 

KG/SEC 

0.0514 

0.0514 

0.0 

0.0 

0.0 

FURFURAL 

KG/SEC 

0.0353 

0.0353 

0.0 

0.0 

0.0 

TOTAL 

KG/SEC 

77.6527 

0.0868 

77.5658 

45.4700 

38.9500 

TEMP 

K 

306. 1500 

306. 1500 

306. 1500 

300. 1500 

826.9415 

PRES 

N/SQM 

. 10130+06 

. 10130+06 

. 10130+06 

. 10130+06 

. 14000+07 

ENTHALPY 

J/KG 

15919+00- 

.23134+07- 

•. 1 3933+08- . 13962+08- 

. 12391+08 

VFRAC 


0.0 

0.0 

0.0 

0.0 

1 . 0000 

LFRAC 


1 . 0000 

1 . 0000 

1 . 0000 

1 . 0000 

0.0 

ENTROPY 

J/KG-K 

-9170.5223- 

3490 . 7728-9 178.331 4-9267 . 7333- 

1674.5832 

DENSITY 

KQ/CUM 

843.4797 

1162.9854 

843.2512 

846.6394 

3.7048 

AV6 MW 


18.0319 

111. 8070 

18.0150 

18.0150 

18.0150 


SUBSTREAM! NCPSD 
CCLOSE KG/SEC 

structure: 

0.0 

NON CONVENTIONAL 
0.0 0.0 

19.6400 

0.0 

ACLOSE 

KQ/SEC 

0.0 

0.0 

0.0 

7.2400 

0.0 

CLBOSE 

KG/SEC 

0.0 

0.0 

0.0 

0.0 

0.0 

GCOSE 

KG/SEC 

8.6022 

0.0 

8.6022 

0.0 

0.0 

ENZYME 

KQ/SEC 

0.0 

0.0 

0.0 

0.0 

0.0 

HCLOSE 

KG/SEC 

0.0 

0.0 

0.0 

10.7100 

0.0 

XYLOSE 

KG/SEC 

9.9989 

0.0 

9.9989 

0.0 

0.0 

LIGNIN 

KG/SEC 

0.0 

0.0 

0.0 

7.8800 

0.0 

TOTAL 

KG/SEC 

18.6012 

0.0 

18.6012 

45.4700 

0.0 

TEMP 

K 

306. 1500 

MISSING 

306. 1500 

300. 1500 

MISSING 

PRES 

N/SOM 

.10130+06 . 

10130+06 

. 10130+06 

. 10130+06 

MISSING 

ENTHALPY 

J/KQ 

. 70460+06 

MISSING 

. 70460+06 

. 69079+06 

MISSING 

VFRAC 


0.0 

MISSING 

0.0 

0.0 

MISSING 

LFRAC 


0.0 

MISSING 

0.0 

0.0 

MISSING 

DENSITY 

KG/CUM 

705.4500 

MISSING 

705.4500 

705.4500 

MISSING 

AVQ MW 


1 . 0000 

1 . 0000 

1 . 0000 

1 . 0000 

1.0000 

COMPONENT ATTRIBUTES! 
ENZYME CAUSRl 

ELEMl 

MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM2 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM3 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM4 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEMS 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM6 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM7 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEMS 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEM9 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

ELEMIO 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 


SUBSTREAM ATTRIBUTES! 
PSD 
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THE ABE PROCESS 
STREAM SECTION 


15 16 17 1 4 (CONTINUED) 


STREAM ID 

15 

16 

17 

1 

4 

FRACl 

0. 3289 

0.3289 

0.3289 

0.0 

MISSING 

FRAC2 

0.1451 

0.1451 

0.1451 

0.0 

MISSING 

FRAC3 

0. 1266 

0. 1266 

0.1266 

0.0 

MISSING 

FRAC4 

0. 1205 

0. 1205 

0. 1205 

0.0 

MISSING 

FRAC5 

0.2542 

0.2542 

0.2542 

0.0 

MISSING 

FRAC6 

0.0244 

0.0244 

0.0244 

0.0 

MISSING 

FRAC7 

0.0 

0.0 

0.0 

0.0 

MISSING 

FRAC0 

0.0 

0.0 

0.0 

0.0 

MISSING 

FRAC9 

0.0 

0.0 

0.0 

0. 5000 

MISSING 

FRACIO 

0.0 

0.0 

0.0 

0.5000 

MISSING 
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THE ABE PROCESS 
STREAM SECTION 

7 

STREAM ID 7 

FROM : 

TO : ACIDIN 

class: mixncpsd 


STRUCTURE: CONVENTIONAL 
0.0 
1 . 9SB7 
0.0 
0.0 
1 . 9SB7 
7S9.37S0 
. 14000+07 
. B0399+06 
1 . 0000 
0.0 

158. 1272 
21.7546 
98.0800 


substream; mixed 

WATER KQ/SEC 
H2S04 KQ/SEC 
HMF KQ/SEC 

FURFURAL KQ/SEC 
TOTAL KQ/SEC 
TEMP K 
PRES N/SQM 
ENTHALPY J/KQ 
VFRAC 
LFRAC 

ENTROPY J/KQ-K 
DENSITY KQ/CUM 
AVG MW 


K 

N/SQM 


J/KQ-K 

KQ/CUM 


substream: ncpsd 

CCLOSE KQ/SEC 
ACLOSE KG/SEC 
CLBOSE KQ/SEC 
QCOSE KQ/SEC 
ENZYME KQ/SEC 
HCLOSE KG/SEC 
XYLOSE KQ/SEC 
LIQNIN KQ/SEC 
TOTAL KQ/SEC 
TEMP K 

PRES N/SQM 
ENTHALPY J/KQ 
VFRAC 
LFRAC 
AVQ MW 

COMPONENT ATTRIBUTES: 

ENZYME CAUSRl 

ELEMl 

ELEM2 

ELEM3 

ELEM4 

ELEM5 

ELEM6 

ELEM7 

ELEM8 

ELEM9 

ELEMIO 

SUBSTREAM ATTRIBUTES: 

PSD 

FRACl 


structure: non conventional 
0.0 
0.0 
0.0 
0-0 
0.0 
0.0 
0.0 
0.0 
0.0 

MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

1.0000 


MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

MISSING 

MISSING 


MISSING 
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THE ABE PROCESS 
STREAM SECTION 


QHTR QACHYO QQNCH QCOOLR 

STREAM ID QHTR QACHYD QQNCH QCOOLR 

FROM : HEATER 1 ACIDHYD QUENCH COOLER 

TO : 

CLASS! HEAT HEAT HEAT HEAT 

STREAM attributes: 

HEAT 

Q WATT -.11547+09-. 15400+08 .20338+09 .30551+08 


105 



M.E.T. 


SHAFT 1 

STREAM 
FROM : 
TO : 

class: 

STREAM 

WORK 

P 


/ VAX ASPEN SYSTEM RELEASE 6 ♦♦ PROD VERSION ♦♦ DATE: U/11/83 PAGE 37 

THE ABE PROCESS 
STREAM SECTION 


ID SHAFTl 

PUMPl 

WORK 

ATTRIBUTES: 

WATT. - 10731+06 
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